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It  is  observed  that  topography  plays  an  important  role  in  forcing  planetary  scale 
stationary  waves  as  well  as  influencing  synoptic  scale  disturbances.  The  purpose  of  this 
study  is  to  investigate  topographic  influences  on  planetary  and  synoptic  scale  waves  in  a 
baroclinic  fluid  using  the  quasi-geostrophic  potential  vorticity  equation  with  wavenumber 
two  bottom  topography  and  a  steady,  wavy  basic  state  on  an  f-plane. 

The  steady  basic  state  is  sinusoidal  in  x  with  linear  vertical  shear  and  represents  an 
exact  solution  to  the  quasi-geostrophic  potential  vorticity  equation  and  thermodynamic 
equations  at  the  upper  and  lower  boundaries.  The  amplitude  of  the  wavy  basic  is  highly 
dependent  on  forcing  by  a  constant  zonal  flow,  UD.  UD  is  a  parameter  used  to  represent 
momentum  flux  convergence  into  the  mid-latitudes. 

The  results  from  a  linear  analytical  solution  show  that  any  mode  with  wavenumber  n 
which  interacts  with  the  two-wave  topography  will  introduce  two  additional  modes,  n±2. 
All  three  modes  will  have  the  same  frequency.  The  results  of  a  linear  stability  analysis 
from  a  numerical  model  show  that  the  growth  rate  for  the  most  unstable  mode  is  not  highly 
dependent  on  the  mountain  height  or  zonal  flow,  UQ.  When  there  is  a  wavy  basic  state 
forced  by  U0,  the  perturbations  are  forced  to  follow  the  basic  state.  When  UD  is  zero, 
there  is  no  wavy  basic  state;  however,  the  perturbations  still  are  forced  to  wave  over  the 
mountain. 
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CHAPTER  1 


INTRODUCTION 

During  the  past  decade,  there  has  been  considerable  research  into  the  role  of 
topography  on  planetary  and  synoptic  scale  flow.  The  research  has  been  observational, 
experimental  and  theoretical.  In  this  thesis,  we  review  some  observational  and 
experimental  research  and  present  a  simple  theoretical  model  of  baroclinic  flow  over 
topography.  We  look  to  the  observational  evidence  for  guidance  in  developing  the 
theoretical  model.  Wavenumber  two  dominates  the  nearly  stationary  pattern  in  the 
atmosphere  and  the  annulus.  So,  we  chose  a  steady  wavy  basic  state  forced  by 
wavenumber  two  topography.  On  the  smaller  scale,  baroclinic  disturbances  tend  to  follow 
the  storm  tracks  as  the  flow  over  the  Rocky  Mountains.  Therefore,  we  did  a  linear  stability 
analysis  of  the  wavenumber  two  basic  state  to  investigate  the  relationship  between 
disturbances  and  the  basic  state. 

Blackmon  (1976),  Blackmon  et  al.  (1977),  Blackmon  etal.  (1979),  and  Wallace 
(1983),  Blackmon  etal.  (1984),  and  Grose  and  Hoskins  (1979)  studied  the  Northern 
Hemisphere  circulation  and  focused  on  traveling  waves  of  various  spatial  and  temporal 
scales.  From  their  studies  of  the  500  mb  geopotential  height  field,  they  found  wavenumber 
two  dominates  and  is  nearly  stationary  with  high  pressure  ridges  over  the  major  mountain 
chains.  Lau  (1979),  Held  (1983),  and  Wallace  (1983)  studied  the  climatology  of  the 
Northern  Hemisphere.  Again,  the  vertical  and  horizontal  cross  sections  for  the  winter 
months  showed  high  pressure  ridges  over  the  major  terrain  features  with  wavenumber  two 
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dominating.  Grose  and  Hoskins  (1979)  showed  that  wavenumber  two  dominates  the 
Earth’s  topographic  features.  Altogether,  these  studies  show  that  the  Northern 
Hemisphere  large  scale  circulation  is  strongly  influenced  by  the  Earth’s  orographical 
features. 

Observational  studies  of  synoptic  scale  transient  eddies  and  the  storm  tracks  in  the 
Northern  Hemisphere  show  that  they  are  influenced  by  topography.  Blackmon  (1976), 
Blackmon  etal.  (1977),  Wallace  (1983),  Blackmon  etal.  (1979),  Hsu  (1987),  and  Buzzi 
and  Tosi  (1989)  completed  statistical  analyses  of  high  frequency  traveling  eddies.  Hsu  and 
Buzzi  and  Speranza  showed  that  traveling  eddies  tend  to  move  north  and  rotate 
anticyclonically  as  they  flow  over  the  Rocky  Mountains.  By  using  the  variance  of  the  data, 
they  showed  that  the  storm  track  is  interrupted  by  the  mountain.  There  is  strong  evidence 
of  the  storm  track  over  the  Pacific  Ocean  and  east  of  the  Rocky  Mountains  but  there  is  a 
minimum  in  activity  over  the  mountains. 

Pfeffer  etal.  (1989)  and  Li  etal.  (1986)  conducted  experimental  studies  using  a 
thermally  driven,  rotating  annulus  with  two  wave  bottom  topography  which  simulates  the 
behavior  of  the  large  scale,  atmospheric  flow.  They  also  found  there  is  a  standing 
wavenumber  two,  with  high  pressure  associated  with  the  mountain  ridge.  They  showed 
that  as  the  rotation  rate  increases,  the  high  pressure  ridge  of  the  standing  wave  is  shifted 
upstream  from  the  mountain. 

Some  theoretical  work  has  been  done  using  wavy  basic  states  which  may  be  forced  by 
topography  or  other  atmospheric  phenomena.  Niehaus  (1980)  used  a  wavy  basic  state  that 
looked  like  the  observed  climatology  for  winter  to  do  a  linear  stability  analysis  of  a 
continuous,  baroclinic  model.  She  compared  wavy  basic  states  and  skewed  basic  states 
which  have  zonal  as  well  as  meridional  structure.  She  found  the  skewed  basic  state  to  be 
more  favorably  configured  for  unstable  waves  to  grow.  Derome  (1984)  showed  how  to 
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find  a  steady-state,  wavy  basic  state  for  a  baroclinic  model  that  is  an  exact  solution  to  the 
nonlinear  quasi-geostrophic  potential  vorticity  equation  and  thermodynamic  conditions  at 
the  upper  and  lower  boundaries.  The  wavy  basic  state  indicates  the  storm  track  regions 
and  regions  of  localized  growth. 

A  considerable  amount  of  theoretical  work  has  been  done  on  the  subject  of  topographic 
instability.  Topographic  instability  may  be  defined  as  a  process  of  wave  growth  without 
propagation.  It  is  caused  by  the  interaction  of  the  mean  zonal  flow  with  topography.  The 
form  drag  equation  couples  the  acceleration  of  the  mean  zonal  flow  with  the  drag  induced 
by  the  topography.  Chamey  and  DeVore  (1979)  used  a  barotropic  model  and  Chamey  and 
Straus  (1980)  used  a  baroclinic  model  to  investigate  the  phenomenon  of  topographic 
instability.  Revell  and  Hoskins  (1984)  argued  that  topographic  instability  can  occur  in 
Chamey  and  DeVore’s  barotropic  model  but  probably  does  not  exist  in  Chamey  and 
Straus’  baroclinic  model  since  more  realistic  primitive  equation  models  do  not  show  this 
form  of  instability.  They  suggest  that  the  lack  of  modes  in  the  meridional  direction  in 
Chamey  and  Straus’  model  was  the  cause.  Kallen  (1983)  showed  that  topographic 
instability  can  occur  in  a  baroclinic  model  like  Chamey  and  Straus’  if  a  source  of 
momentum  is  included  to  drive  a  zonal  flow.  Buzzi  et  al.  (1984)  also  used  a  baroclinic 
model  with  a  barotropic  forcing  component  to  show  topographic  instability  of  the  type 
found  by  Chamey  and  DeVbre. 

In  this  thesis,  we  developed  a  simple  mathematical  model  to  describe  a  planetary  scale 
basic  flow  and  the  smaller  scale,  baroclinic  disturbances.  We  scaled  the  governing 
equations  so  that  they  would  describe  these  two  different  scales  of  motion.  Based  on  the 
scaling  of  the  governing  equations,  we  found  that  the  flows  could  be  described  using  the 
quasi-geostrophic  potential  vorticity  equation  and  thermodynamic  equations  at  the  upper 
and  lower  boundaries.  The  effects  of  topography  enter  the  problem  through  the 


thermodynamic  equation  at  the  lower  boundary.  The  equations  are  valid  for  non- 
dimensional  mountain  heights  on  the  order  of  the  Rossby  number. 
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We  linearized  the  model  about  a  basic  state  that  was  wavenumber  two  in  the  x-direction 
and  had  linear  vertical  shear.  We  took  the  model  to  be  periodic  in  x  and  y  with  a  rigid 
upper  lid  and  cos(2Kx)  bottom  topography.  Since  the  model  is  periodic  in  y,  side  wall 
boundary  conditions  are  not  considered. 

To  force  the  x-dependent  basic  state,  we  need  a  constant  zonal  flow  at  the  surface 
which  interacts  with  the  topography.  In  the  atmosphere,  a  constant  zonal  flow  results  from 
eddy  momentum  flux  convergence  into  the  mid-latitudes.  Since  we  chose  periodic  sidewall 
boundary  conditions,  the  model  cannot  describe  momentum  flux  convergence.  As  a  result, 
the  constant  zonal  flow  enters  this  problem  as  a  free  parameter. 

We  solved  the  linear  quasi-geostrophic  potential  vorticity  equation  and  thermodynamic 
equations  at  the  upper  and  lower  boundaries  by  analytical  and  numerical  techniques.  The 
analytical  solution  uses  asymptotic  expansions  in  a  small  parameter  related  to  the  mountain 
height.  The  leading  order  problem  is  the  Eady  problem.  At  the  next  order  enters  the 
correction  to  the  Eady  problem  due  to  topography  and  wavy  basic  state.  We  find  that  the 
frequency  from  the  leading  order  problem  is  also  the  frequency  for  the  next  order  where  the 
correction  due  to  topography  enters  the  problem.  Thus,  all  of  the  modes  found  in  the  first 
two  orders  of  this  problem  will  have  the  same  frequency  and  growth  rate.  The  numerical 
solution  involves  Fourier  transforming  the  equation  from  grid  space  to  spectral  space.  This 
gives  a  set  of  amplitude  equations  which  can  be  solved  using  an  eigenvector,  eigenvalue 
routine  to  find  growth  rates,  frequencies,  and  amplitudes  of  the  spectral  coefficients.  The 
eigenvectors  are  used  to  reconstruct  the  streamfunction.  The  results  from  the  analytical 
model  help  us  to  understand  the  results  from  the  numerical  model. 

In  Chapter  2,  we  review  observational  evidence  from  the  atmosphere  and  the  annulus 
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that  suggests  the  influence  of  topography  on  the  flow.  In  Chapter  3,  we  develop  the 
baroclinic  model  we  wili  use  to  investigate  flow  over  topography.  In  Chapter  4,  we  solve 
the  model  using  a  numerical  and  an  analytical  technique.  Finally,  in  Chapter  5,  we  discuss 
results  and  conclusions. 


CHAPTER  2 

ATMOSPHERIC  AND  ANNULUS  OBSERVATIONS 

Observations  from  the  atmosphere  and  from  annulus  experiments  with  topography 
show  flow  patterns  that  are  related  to  topographic  forcing.  By  studying  atmospheric 
variability  we  may  learn  about  the  geographical  distribution  of  the  geopotential  height  field 
at  different  vertical  levels  for  different  time  scales.  Low  frequency  variability  indicates  the 
nearly  stationary  pattern  in  the  atmosphere.  The  long  wave  spectrum  of  low  frequency 
variability  dominates  the  climatology  of  the  500  mb  geopotential  height  field.  From  the 
climatology  we  find  the  relationship  between  the  nearly  stationary,  planetary  scale  waves 
and  topographic  features.  High  frequency  variability  indicates  the  storm  tracks  and  regions 
of  high  eddy  activity.  Thus,  we  can  learn  about  the  relationship  of  high  frequency  eddies 
to  topographic  features  and  to  the  long  wave  environment.  Laboratory  experiments  using  a 
rotating  annulus  with  a  thermally  driven  fluid  and  topography  provide  another  way  to 
understand  the  interaction  of  large  scale  flow  with  topography.  In  this  chapter,  we  look  at 
observational  evidence  of  topographic  influences  on  atmospheric  and  annulus  flows. 


^wiTTmoTn  iTTifmt J  IwTITnrnTiT 


If  topography  plays  an  important  role  in  atmospheric  flows,  then  it  should  leave  a 
signature  that  can  be  identified  by  observational  studies  of  atmospheric  variability  of  the 
Northern  Hemisphere.  There  should  be  a  signature  due  to  topography  in  the  nearly 
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stationary,  long  wave  pattern  as  well  as  the  higher  frequency  baroclinic  eddies.  Blackmon 
(1976)  and  Blackmon  etal.  (1977,  hereafter  BWLM)  studied  the  geographical  distribution 
of  low  frequency  variability  on  the  time  scale  of  greater  than  10  days.  Low  frequency 
variability  indicates  the  nearly  stationary  flow  pattern.  They  analyzed  data  from  the  500  mb 
NMC  analysis  for  18  winter  seasons  in  the  Northern  Hemisphere.  Based  on  the  time 
filtering  scheme  they  chose,  a  low  pass  filter  gathered  data  with  periods  greater  than  10 
days,  a  medium  pass  or  band  pass  filter  gathered  data  with  periods  between  2.5  -  6  days. 

Figure  2.1  from  Blackmon  (1976)  shows  maps  of  the  500  mb  rms  geopotential  height 
surface  for  the  low  pass  filtered  data  for  all  wavelengths.  The  pattern  is  almost  identical  to 
the  unfiltered  500  mb  rms  geopotential  surface  (not  shown).  Thus,  the  500  mb 
climatology  shown  in  Figure  2.2  from  Blackmon  (1976)  is  dominated  by  low-frequency 
components.  The  500  mb  climatology  shows  a  high  pressure  ridge  over  the  Rocky 
Mountains  and  the  Tibetan  Plateau.  Wavenumber  two  dominates  the  pattern.  In  a 
topographical  study  of  the  Earth’s  surface,  Grose  and  Hoskins  (1979)  used  a  spectral 
analysis  to  show  that  wavenumber  two  dominates  the  surface. 

There  are  many  possible  mechanisms  in  addition  to  topography  which  could  force  low 
frequency  variability.  One  explanation  suggests  it  is  due  to  slowly  varying  thermal  forcing 
from  the  oceans.  Another  explanation  uses  the  index  cycle  in  which  there  are  vacillations 
from  the  high  index  state  with  strong  westerlies  and  weak  eddies  to  the  low  index  state 
with  weak  westerlies  an  large  amplitude  waves.  A  third  explanation  for  low  frequency 
variability  is  the  interaction  of  orographically  forced  waves  with  slow  moving  transient 
planetary  waves  such  as  Rossby  waves  which  interact  and  produce  long  lived,  large 
amplitude  waves.  A  fourth  explanation  is  forcing  by  high  frequency  baroclinic  eddies 
which  transfer  their  energy  up  the  scale  to  the  longer  waves. 

Finally,  a  fifth  explanation  suggests  that  topography  forces  the  low  frequency 
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Figure  2.1  Map  of  low  pass  filtered  rms  geopotential  at  500  mb  for  nine  winters. 
Emphasizes  fluctuations  in  periods  longer  than  10  days.  The  contribution  from 
the  mean  annual  cycle  has  been  removed.  Contour  interval  is  10  m.  From 
Blackmon  (1976). 


Figure  2.2  Climatological  mean  500  mb  height  field  for  winter.  Contour  interval  is  50 
m.  From  Blackmon  (1976). 
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variability  of  the  atmosphere.  Figure  2.3  taken  from  Blackmon,  Lee,  and  Wallace  (1984, 
hereafter  BLW)  shows  (a)  the  500  mb  and  (b)  the  1000  mb  height  patterns  averaged  over 
24  days  beginning  on  31  December  1980.  For  this  24  day  period,  a  strong  ridge  is 
centered  over  the  Rocky  Mountains  with  deep  troughs  over  the  Pacific  and  North  America. 
This  pattern  has  been  observed  during  many  winter  periods.  When  the  remaining  days  of 
that  winter  were  averaged  in,  the  composite  map  looked  very  much  like  the  500  mb 
climatological  map  shown  in  Figure  2.2.  Since  this  pattern  has  been  observed  many  times, 
it  may  be  said  that  this  represents  one  of  the  atmosphere’s  preferred  responses  to  the  Rocky 
Mountains. 

However,  there  are  other  ways  the  atmosphere  could  adapt  to  the  Rocky  Mountains. 
Figure  2.4  also  taken  from  BLW  shows  (a)  the  500  mb  and  (b)  the  1000  mb  geopotential 
height  patterns  for  22  days  during  the  winter  of  1980.  In  this  figure,  the  500  mb  ridge  is 
located  over  the  Gulf  of  Alaska  so  that  the  flow  runs  parallel  to  the  Rockies.  Figure  2.3 
showed  the  flow  was  nearly  perpendicular  to  the  Rocky  Mountains.  At  1000  mb,  the  flow 
shown  in  Figure  2.4  is  easterly  over  the  Rockies  whereas  it  is  westerly  in  Figure  2.3. 

Also,  Figure  2.3  shows  the  vertical  structure  is  more  baroclinic  than  in  Figure  2.3  where 
the  wave  is  nearly  vertically  stacked.  In  Section  3.6,  we  will  consider  a  baroclinic  basic 
state  which  is  torced  by  topography.  For  this  basic  state,  the  flow  will  be  similar  to  the 
case  in  Figure  2.3. 

2.1.2  High  Frequency  Variability 

Blackmon  (1976),  BWLM,  Wallace  and  Blackmon  (1983),  Hsu  (1987),  and  Buzzi  and 
Tosi  (1989)  studied  the  geographical  distribution  of  high  frequency  variability.  High 
frequency  variability  indicates  storm  track  regions.  Hsu  (1987)  applied  statistical  analysis 
on  band  pass  filtered  data  in  the  vicinity  of  mountain  ranges.  He  found  that  synoptic  scale 
eddies,  both  cyclones  and  anticyclones,  in  the  vicinity  of  large  mountain  ranges  tend  to 
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Figure  2.3  Composite  geopotential  height  patterns  for  24-day  period  beginning  31 
December  1980,  based  upon  daily  analyses  made  at  ECMWF.  (a)  500  mb  height,  contour 
interval  is  60  m  and  (b)  10  mb  height,  contour  interval  30  m.  From  Blackmon,  Lee,  and 
Wallace  (1984). 


Figure  2.4  Composite  geopotential  height  patterns  for  the  22  days  1-9  and  18-24 
December  1980  and  5-10  February  1981,  inclusive.  Based  on  daily  analyses  made  at 
ECMWF.  From  Blackmon,  Lee,  and  Wallace  (1984). 
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rotate  anticyclonically  around  the  mountain.  This  may  be  understood  in  terms  of  a 
barotropic  explanation  which  says  that  potential  vorticity  is  conserved.  The  potential 
vorticity  may  be  expressed  as 

o  _  C +  fo  +  Poy 

V  D  +  h 

where  Q  is  the  potential  vorticity,  fQ  is  the  planetary  vorticity  evaluated  at  some  mid¬ 
latitude,  P0  is  the  beta-effect  at  that  latitude,  h  is  the  mountain  height,  D  is  the  average 
depth  of  the  fluid,  and  £  is  the  relative  vorticity.  As  the  eddy  flows  up  the  mountain  slope, 
the  depth  of  the  fluid  decreases.  To  conserve  absolute  vorticity,  the  relative  and/or 
planetary  vorticity  also  must  decrease.  Since  it  is  observed  that  the  eddies  move  to  the 
north  as  they  flow  over  the  mountains,  the  planetary  vorticity  also  increases.  Therefore, 
the  relative  vorticity  must  decrease  and  the  eddy  rotates  anticyclonically.  A  baroclinic 
explanation  for  the  eddies  turning  north  would  be  in  terms  of  adiabatic  warming  and 
cooling  of  the  fluid.  The  adiabatic,  thermodynamic  energy  equation  could  be  used  to  show 
that  as  the  fluid  rises,  it  cools  adiabatically.  The  adiabatic  cooling  would  induce  a 
meridional  circulation  which  advects  warm  air  from  the  south  and  also  advects  the  eddy  to 
the  north.  In  the  atmosphere,  there  must  be  a  combination  of  both  the  baroclinic  and 
barotropic  mechanisms.  In  Chapter  4  of  this  thesis,  we  will  use  a  baroclinic  model  with  a 
constant  zonal  flow  to  show  the  northward  excursion  of  eddies  as  they  pass  over  a 
mountain;  in  particular  we  propose  to  consider  the  f-plane  dynamics. 

Buzzi  and  Tosi  (1989,  hereafter  BT)  extended  the  work  of  Hsu.  They  examined  the 
spatial  structure  of  high-frequency  variability  near  mountains  including  the  Rocky 
Mountains.  In  particular,  they  studied  the  effect  of  orography  on  baroclinic  eddies  that 
have  periods  of  a  week  or  less.  They  used  time  series  of  geopotential  height  at  850  and 
500  mb  derived  from  6-hour  initialized  analyses  of  ECMWF  for  the  period  November 
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1979  -  December  1986.  They  analyzed  the  time  series  with  a  band  pass  filter  to  isolate  the 
baroclinic  eddies  which  have  periods  less  than  one  week.  They  included  data  from  all 
seasons  but  high-frequency  variability  is  greatest  during  winter  so  their  results  should  be 
more  representative  of  the  colder  months. 

Figure  2.5  from  BT  shows  a  map  of  (a)  the  850  mb  mean  geopotential  height,  (b)  the 
850  mb  standard  deviation  of  the  high-pass  filtered  geopotential,  (c)  the  500  mb  mean 
geopotential  height,  and  (d)  the  500  mb  standard  deviation  of  the  high-pass  filtered 
geopotential  over  the  Rocky  Mountains.  The  relative  maximum  in  the  standard  deviation 
over  the  Pacific  to  the  west  of  the  mountains  shows  the  storm  tracks  at  both  850  and  500 
mb.  The  storm  track  is  interrupted  over  the  mountains  where  the  standard  deviation  shows 
a  minimum  at  both  levels.  To  the  east  of  the  mountains,  there  is  a  relative  maximum  at  850 
mb  that  suggests  the  storm  track  migrates  to  the  south  at  low  levels  but  simply  reforms  east 
of  the  mountains  at  500  mb. 

In  addition  to  studying  the  variance  of  the  eddies,  BT  used  an  analysis  technique  based 
on  one-point  lag  correlation  statistics.  Their  technique  was  to  select  a  base  gridpoint  near 
the  mountain  range.  Then,  they  computed  lag-correlations  between  the  time  series  at  the 
base  gridpoint  and  the  time  series  at  all  other  gridpoints  and  at  all  levels.  The  lag 
correlation  maps  provide  information  on  the  spatial  structure  of  the  eddies  in  the  region 
surrounding  the  base  gridpoints.  This  is  possible  since  eddies  preserve  spatial  structure 
over  a  time  period  greater  than  its  local  period  and  over  a  distance  greater  than  its  mean 
wavelength.  Lag  correlation  maps  may  be  compared  to  synoptic  maps  but  heights  cannot 
be  interpreted  as  the  geopotential  height. 

Figure  2.6  from  BT  shows  the  one-point  lag  correlation  maps  for  several  lag  times  of 
the  high-pass  filtered  data  of  the  sea-level  pressure  at  (54°N,129°W)  marked  with  a  dot 
and  the  sea-level  pressure  at  all  other  gridpoints.  To  produce  the  maps,  they  selected  data 


Figure  2.5  Climatology  of  the  850  and  500  mb  levels  for  the  area  of  the  Rocky 
Mountains,  (a)  850  mb  mean  geopotential  height  and  (b)  standard  deviation  of  the  high 
pass  filtered  geopotential  height  divided  by  the  sine  of  the  latitude,  (c)  500  mb  mean 
geopotential  and  (d)  is  the  standard  deviation.  From  Buzzi  and  Tosi  (1989). 
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Figure  2.6  Maps  of  lagged  correlation  coefficients  between  the  high-pass  filtered  850mb 
level  at  54°N,  129°W  and  the  high  pass  filtered  850  mb  level  at  every  other  grid  point. 

The  black  dot  show  this  location.  Each  plot  shows  a  different  lag  time,  (a)  -24  hr,  (b)  -12 
hr,  (c)  zero  lag,  (d)  +12  hr,  (e)  +24  hr,  (f)  +36  hr.  Contour  interval  is  0.05  and  negative 
values  are  dashed.  From  Buzzi  and  Tosi  (1989). 


Figure  2.7  Correlation  values  between  high-pass  filtered  850  mb  level  and  high  pass 
filtered  500  mb  level  at  every  grid  point.  Area  includes  the  Rocky  Mountains.  From  Buzzi 
and  Tosi  (1989). 
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only  if  the  wind  direction  was  within  an  angle  of  ±10°  of  the  zonal  wind  direction.  Figure 
2.6  shows  that  the  path  followed  by  the  eddies  is  anticyclonically  curved  around  the 
mountain  as  noted  by  Hsu  (1987). 

Figure  2.7  from  BT  shows  the  correlation  between  the  850  mb  and  500  mb  levels.  The 
two  levels  are  highly  correlated  to  the  west  of  the  Rocky  Mountains  but  there  is  a  minimum 
to  the  east  of  the  mountains.  This  may  be  interpreted  as  the  local  maximum  tilt  of  the 
vertical  eddy  structure.  The  maximum  vertical  tilt  is  due  to  low  level  temperature 
anomalies;  as  the  air  flows  down  the  mountains  slope,  it  warms  adiabatically. 

IAmL  -Climatology  of  Northern  Hemisphere  Wintertime  Stationary  Waves 

Held  (1983)  ran  a  general  circulation  model  with  and  without  topography  and  showed 
that  topography  plays  a  major  role  in  determining  the  longitude  of  the  major  wintertime 
troughs  and  ridges  in  the  Northern  Hemisphere  500  mb  geopotential  height  pattern.  The 
climatology  for  the  Northern  Hemisphere  winter  shows  a  similar  longitude  dependence. 
This  suggests  that  topography  plays  an  important  role  in  determining  the  Northern 
Hemisphere  climatology. 

In  a  climatological  study,  Lau  (1979)  showed  both  horizontal  and  vertical  cross 
sections  of  geopotential  that  appear  to  be  strongly  influenced  by  the  topography.  Figure 
2.8  (a)  -  (c)  from  Lau  shows  a  longitude-height  section  of  the  stationary  wave  geopotential 
height  pattern  at  three  different  latitudes  for  the  winter  season,  derived  from  1 1  years  of 
NMC  operational  analyses.  A  two  wave  pattern  dominates  the  cross  section  at  45°N  and 
60°N.  There  are  high  pressure  ridges  over  the  Rocky  Mountains  (~1 10°W)  and  the 
Tibetan  Plateau  (~115°E).  He  showed  that  these  standing  features  appear  in  individual 
years  as  well  as  the  climatology. 

Figure  2.2  shows  the  climatological  mean  for  the  500  mb  geopotential  height  field  for 
January.  Again,  wavenumber  two  dominates  the  geopotential  height  pattern  and  ridges  are 
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located  over  the  Rocky  Mountains  and  the  Tibetan  Plateau.  The  sea-level  pressure  field 
(not  shown)  does  not  show  a  wavenumber  two  forced  by  mountains  but  rather  strong 
thermal  influences  along  the  coast  Land-sea  temperature  differences  play  a  major  role  in 
the  low  level  climatology.  Held’s  (1990)  numerical  model  also  produces  this  observed 
climatological  pattern. 


2u2 _ Annulus  Observations 

To  simulate  some  of  the  processes  of  the  general  circulation  of  the  earth’s  atmosphere, 
Pfeffer  etal.  (1989)  and  Li  etal.  (1986)  conducted  numerous  laboratory  experiments  in  a 
thermally  driven,  rotating  annulus  with  flow  forced  over  a  two-wave  sinusoidal  bottom 
topography.  The  selection  of  wavenumber  two  is  motivated  by  the  dominance  of  the 
wavenumber  two  in  the  earth’s  topography.  Results  from  annulus  experiments  may  be 
compared  to  theoretical  models  intended  to  describe  flow  in  the  annulus  over  a  broad  range 
of  parameters.  In  Chapter  4  of  this  thesis,  we  present  more  details  about  the  experiment 
and  the  parameter  ranges  studied. 

Figure  2.9,  taken  from  Pfeffer  et  al.  (1989),  shows  the  zonal  wavenumber  spectra  of 
streamfunction  at  mid-radius  in  the  annulus  for  several  rotation  rates.  Figure  2.9  (a) 
represents  the  time  average  of  the  sequence  of  individual  spectra  of  streamfunction  which 
includes  traveling  and  standing  waves.  At  low  rotation  rates,  wavenumbers  two  and  five 
dominate  but  as  the  rotation  rate  increases,  wavenumber  four  dominates  and  the  even 
harmonics  have  the  largest  amplitudes.  Figure  2.9  (b)  represents  the  spectra  of  the  time 
mean  of  the  streamfunction.  This  figure  shows  the  climatology  or  the  standing  waves. 

The  amplitude  of  wavenumber  two  dominates  the  climatology  and  the  even  harmonics  have 
the  largest  amplitudes.  The  fact  that  the  even 
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Figure  2.9  Zonal  wavenumber  spectra  of  y  at  mid-radius  for  seven  different  rotation 
rates,  (a)  Time  averages  of  the  sequence  of  individual  spectra  of  y  and  (b)  spectra  of  the 
time-averaged  (or  climatological)  y  field  representing  the  standing  waves.  From  Pfeffer 
etal.  (1989). 
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harmonics  dominate  suggest  there  are  non-linear  interactions  with  wavenumber  two  forced 
by  topography.  Figure  2.10,  also  taken  from  Pfeffer  et  al.  (1989),  shows  the 
climatological  time-averaged  and  zonally-averaged  velocity  profile,  U,  as  a  function  of 
radius,  the  time-averaged  streamfunction,  and  the  time-averaged  or  standing  wavenumber 
two  at  various  rotation  rates.  At  low  rotation  rate,  the  zonal  current  is  strong  and  the 
standing  wave  is  weak.  At  higher  rotation  rates,  the  standing  wave  is  stronger  and  the 
zonal  flow  decreases.  At  low  rotation  rates,  the  highs  in  the  wavenumber  two  pattern  are 
nearly  standing  on  the  mountain  ridge  as  is  observed  in  the  atmosphere.  At  higher  rotation 
rates,  the  high  pressure  is  phase  shifted  upstream  from  the  ridge  of  the  mountain. 

In  summary,  topography  plays  an  important  role  in  determining  the  large  scale  flow  as 
well  as  the  synoptic  scale  flow  as  observed  in  the  atmosphere  and  the  annulus.  The 
climatological  or  nearly  stationary  pattern  appears  to  be  dominated  by  a  standing  wave 
number  two  in  both  the  atmosphere  and  the  annulus.  Also,  topography  affects  the 
baroclinic  eddies  as  observed  in  the  atmosphere.  In  the  remainder  of  this  thesis,  we  will 
use  theoretical  models  to  study  the  standing  wavenumber  two  pattern  and  the  effects  of 
topography  and  standing  waves  on  the  traveling  baroclinic  eddies. 


20 


U  f  Tp?  0 

I  .00 

•  7S 

y/Y  .so 
.25 
.00 
1 .00 

.75 

y/Y  .50 

.25 

.00 
I  .00 

.75 

y/Y  .50 

•  2S 

.00 
I  .00 

.75 

y/Y  .S0 

.25 
.00 
I  .00 

.75 

y/Y  .50 

.25 

.00 
I  .00 

.75 

y/Y  .50 

•  2S 

.00 
I  .00 

.75 

y/Y  -50 

.25 
.00 

.0  .1  .2  3 

a  b  c 

Figure  2.10  Velocity  data  based  on  measurements  taken  at  a  depth  of  5.7  cm  in  seven 
experiments  at  different  rotation  rates,  (a)  Time-  and  zonally  -averaged  velocity  profile 
(U).  Here,  y  is  the  distance  from  the  outer  toward  the  inner  wall  and  Y  is  the  g^p  width  of 
the  annulus,  (b)  Time-averaged  streamf unction,  \y.  Contour  interval  is  0.1  cm2  s'1-  (c) 
Time-averaged  wavenumber  2  component  of  the  stream  function,  determined  by  Fourier 
analyses  of  y.  Units  of  U  are  cm  s'  .  Contour  interval  is  0.02  cm2  s'  .  Arrows  indicate 
locations  of  topographic  ridges.  From  Pfeffer  et  al.  (1989). 


CHAPTER  3 


MODEL  FORMULATION 


In  this  chapter,  we  develop  a  mathematical  model  within  the  quasi-geostrophic 
framework  to  describe  the  interaction  of  the  planetary  and  synoptic  scale  baroclinic  flows 
with  a  two-wave  bottom  topography.  We  selected  wavenumber  two  since  it  dominates  the 
Earth’s  topography.  To  gain  insight  into  the  physical  mechanisms  involved  in  the  flow,  we 
use  a  simple  set  of  governing  equations  and  seek  analytical  solutions.  We  transform  the 
equations  to  sigma  coordinates  and  then  to  non-dimensional  form.  After  simplifying  the 
equations  using  order  of  magnitude  arguments,  we  form  the  quasi-geostrophic  potential 
vorticity  equation.  In  our  model,  the  effects  of  topography  enter  the  problem  through  the 
lower  thermodynamic  equation  and  there  is  a  rigid  upper  lid.  Finally,  we  consider  the 
form-drag  equation. 


Pedlosky  (1987)  gives  the  complete  set  of  governing  equations  in  spherical  coordinates 
for  the  atmosphere.  Bird,  Stewart,  and  Lightfoot  (1960)  give  the  governing  equations  in 
cylindrical  coordinates  for  the  annulus.  We  can  simplify  and  write  the  governing  equations 
in  Cartesian  coordinates  if  we  choose  to: 

(1)  Neglect  curvature  effects. 

(2)  Approximate  the  Earth’s  rotational  effect  with  an  f-plane. 

(3)  Assume  a  static  state  exists  in  the  absence  of  motion. 
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(4)  Use  the  Boussinesq  approximation. 

(5)  Make  the  hydrostatic  approximation. 

(6)  Assume  fluid  parcel  motions  are  adiabatic. 

(7)  Use  no  heat  sources  or  sinks. 

(8)  Neglect  Ekman  dissipation. 

(9)  Neglect  internal  dissipation. 

Using  the  list  of  assumptions  given,  we  write  the  governing  equations  for  the  model  as 
follows: 


+  f  kx\?  +  _L  v  P  =  0, 
dt  p0 

(3.1.1) 

ap 

a^  =  -gp’ 

(3.1.2) 

V- V  +  ^  =0, 
dz 

(3.1.3) 

P  =  P  o0  "  a(T  -  T0)), 

(3.1.4) 

and 

e-B 

n 

© 

(3.1.5) 

where 

h  a  a  a  a 

dt  dt  dx  dx  dz 

The  horizontal  momentum  equation  is  given  by  (3.1.1)  where  V  represents  the  horizontal 
velocity  vector.  Equation  (3.1.2)  is  the  hydrostatic  approximation  and  (3.1.3)  is  the 
continuity  equation.  Equation  (3.1.4)  is  the  equation  of  state  for  a  liquid  where  T0  is  a 
reference  temperature  and  a  is  the  coefficient  of  volumetric  expansion;  this  equation 
satisfies  the  Boussinesq  approximation.  Equation  (3.1.5)  is  the  adiabatic  thermodynamic 
energy  equation. 
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In  this  problem  we  consider  flows  of  two  different  length  scales;  waves  forced  by 
topography  which  are  planetary  in  scale  and  growing  baroclinic  eddies  which  are  synoptic 
in  scale.  We  expect  that  the  governing  equations,  (3.1.1)  -  (3.1.5),  can  describe  both  of 
these  scales,  j^et  us  consider  the  scale  for  the  planetary  waves  and  the  validity  of  the 
assumptions  we  applied.  Namely,  how  important  are  the  curvature  of  the  Earth  and  the 
beta  effect  for  planetary  scale  waves?  In  Appendix  A  we  scale  the  horizontal  momentum 
equations  in  spherical  and  cylindrical  coordinates  to  estimate  the  magnitude  of  the  Earth’s 
curvature  compared  to  advective  terms  and  we  consider  the  beta-plane  approximation  to 
estimate  the  magnitude  of  the  beta  effect  compared  to  the  constant  vorticity  taken  in  the  f- 
plane  approximation.  We  summarize  the  results  here.  In  Section  3.3,  we  scale  the 
equations  for  synoptic  scale  motions. 

To  neglect  curvature  effects  in  models  intended  to  describe  planetary  scale  flow  of  the 
atmosphere  or  in  the  annulus,  the  length  scale  of  the  phenomena  should  be  much  smaller 
than  the  Earth’s  radius  or  the  annulus  radius.  If  »o  a  first  approximation  We  consider  the 
Earth’s  surface  to  be  two  waves  in  the  zonal  direction,  then  the  typical  horizontal  length 
scale  of  a  standing  wave  forced  by  the  topography  is  approximately  one  half  the  magnitude 
of  the  Earth’s  radius.  In  the  annulus  experiment,  there  is  a  two  wave  sinusoidal 
topography  which  will  force  a  standing  wavenumber  two.  Again,  the  horizontal  length 
scale  is  approximately  one  half  the  magnitude  of  the  radius  of  the  annulus.  Thus,  curvature 
effects  are  important  to  describe  properly  planetary  scale  flow  forced  by  topography.  By 
including  these  effects,  the  analytical  solution  would  be  very  complicated  and  a  numerical 
solution  would  be  easier.  We  expect  that  the  curvature  terms  will  bring  some  correction  to 
the  simple,  analytical  model  we  solve.  In  our  model,  we  chose  to  neglect  curvature  effects 
in  favor  of  the  advective  terms. 

We  approximated  the  Earth’s  rotational  effect  with  an  f-plane  rather  than  using  a  beta- 
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plane  approximation.  In  the  atmosphere,  the  standing,  planetary  scale  waves  forced  by 
topography  will  have  large  north-south  excursions  and  the  beta  effect  would  be  important. 
In  the  rotating  annulus,  the  bottom  surface  is  not  radially  dependent  but  there  is  a  sloping 
upper  surface  which  simulates  the  beta  effect.  In  either  case,  a  better  approximation  than 
the  f-plane  is  the  beta-plane.  In  Appendix  A,  we  compare  the  f-plane  and  beta-plane 
approximations  and  show  the  magnitude  of  the  beta  effect  is  of  the  same  order  as  the 
Rossby  number  for  planetary  scale  waves.  To  describe  planetary  scale  waves  we  should 
include  the  beta  effect  if  we  include  the  advective  terms  which  are  of  the  same  order.  To 
describe  both  planetary  scale  and  synoptic  scale  motions,  we  need  the  beta  effect  and 
advective  terms.  We  consider  synoptic  scales  in  Section  3.3.  In  this  model,  we  chose  a 
standing  planetary  scale  wave  forced  by  topography  but  not  by  the  beta  effect.  This 
represents  a  serious  deficiency  in  the  model.  However,  the  stability  analysis  in  Chapter  4 
focuses  on  baroclinic  eddies  within  the  long  wave  environment  where  the  beta  effect  is  not 
as  important  for  synoptic  scale  waves  as  for  planetary  scale  waves. 

In  (3.1.2)  we  assume  that  vertical  accelerations  are  small  compared  to  the  rotational 
effects  so  the  hydrostatic  approximation  is  valid.  For  this  approximation  to  remain  valid  in 
the  model  with  mountains,  the  height  of  the  mountain  should  be  less  than  or  equal  to  the 
magnitude  of  the  Rossby  number  (Pedlosky,  1987). 

3.2  Transformation  of  Equations  To  Sigma  Coordinates 

One  way  to  solve  this  problem  would  be  to  assume  that  the  mountain  height  is  small 
compared  to  a  characteristic  vertical  scale.  Then,  we  could  linearize  the  lower  boundary 
condition  and  apply  it  on  the  z  =  0  surface  as  shown  by  Phillips  (1963).  This 
approximation  is  invalid  if  the  mountain  height  is  much  larger  than  the  Rossby  number. 
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Another  possible  way  to  solve  this  problem  would  be  to  transform  the  equations  to  a  terrain 
following  coordinate  system.  Malguzzi  (1987)  showed  that  if  the  topography  is  steep,  they 
had  to  abandon  the  z-coordinate  system  in  favor  of  a  sigma  coordinate  system.  To  solve 
this  problem,  we  transform  our  equations  to  a  terrain  following  or  sigma  coordinate 
system. 

We  define  a  new  vertical  coordinate  such  that 
q(x  z)  =  D(z '  h(x)) 

(D  -  h(x))  f  3.2.11 


where  D  is  the  average  depth  of  the  fluid,  h  is  the  height  of  the  mountain,  z  is  the  height 
above  the  average  mountain  depth,  and  sigma  is  the  new  vertical  coordinate  as  shown  in 
Figure  3.1.  In  this  coordinate  system,  when  z  =  h(x),  a  -  0.0  and  when  z  =  D,  o  =  D. 
For  this  problem,  we  took  the  mountain  to  be  x-dependent  only.  The  Rocky  Mountains 
extend  for  thousands  of  kilometers  in  the  north-south  direction  and  a  good  approximation 
for  the  model  is  to  take  a  mountain  that  is  infinite  in  y.  In  the  annulus,  the  topography  does 
not  have  any  y-structure. 


Figure  3.1  Cross  section  showing  the  relationship  between  the  a-coordinate  and  the  z- 
coordinate.The  shaded  region  indicates  the  topography.  The  surface  a  =  0  is  coincident 
with  the  topography.  The  surface  z  =  0  is  at  the  midpoint  between  the  trough  and  ridge. 
D  is  the  average  fluid  depth. 


To  transform  (3.1.1)  -  (3.1.5)  from  the  Cartesian  coordinate  system  to  the  sigma 
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coordinate  system  defined  above,  we  use  the  following  relations 


1*1 

./*)  +» 

Li90! 

Uxl 

(z)  1  Sx  1  („)  di 

G\  dx  1 

3F  =  9F/aa  \ 
dz  do  \  dx  )  (z) 


(3.2-2) 


(3.2.3) 


Thus,  we  have  a  general  transformation  from  the  z-coordinate  system  to  the  sigma 
coordinate  system.  For  steep  mountains,  the  sigma  surfaces  in  the  lower  levels  would  be 
very  wavy  when  viewed  in  a  z-coordinate  system.  For  the  case  of  steep  mountains,  our 
transformation  as  defined  here  becomes  invalid.  Pielke  (1984)  and  Dutton  (1976)  show  the 
appropriate  transformation  to  a  generalized  non-orthogonal  coordinate  system  in  the  case  of 
steep  mountains.  For  this  problem,  we  used  mountains  with  gentle  slopes. 

Using  the  transformations  (3.2.2)  and  (3.2.3),  we  find 


(3.2.4) 

(3.2.5) 


Also,  we  find  the  following  operators  apply  in  the  sigma-coordinate  system 


SL  = 
dt(a) 


-  -1 


.  d 


-+V-W  +<j^- 


dt 


do 


and 


J{g) 

w  -’M+ISifllrlJs 


(3.2.6) 


(3.2.7) 

Pedlosky  (1987)  argued  that  in  the  absence  of  rotation,  we  expect  there  is  some 
‘standard”  atmosphere  that  is  only  a  function  of  the  vertical  coordinate.  If  we  form  a  global 


average  of  the  pressure,  P,  and  temperature,  T,  and  assume  we  can  write  the  average  as  the 
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static  pressure,  Ps(z),  and  temperature,  Ts(z),  then  we  can  write  the  total  pressure  and 
temperature  as 


P(x,y,z,t)  =  Ps(z)  +  P'(x,y,z,t), 
and 

T(x,y,z,t)  -  T0  =  Ts(z)  +  T(x.y,z,t), 


(3.2.9) 

(3.2.10) 


where  F (x,y,z,t)  and  T (x,y,z,t)  represent  spatial  and  temporal  perturbations  to  the  static 


state  Ps  and  Ts.  When  we  transform  the  equations  to  the  sigma  coordinate  system,  we  will 
keep  Ps(z)  and  Ts(z)  as  functions  of  z. 


Using  these  operators  and  definitions,  we  may  write  the  governing  equations  in  sigma 
coordinates  as  follows 


(3.2.11) 

(3.2.12) 

(3.2.13) 

(3.2.14) 


At  the  upper  and  lower  rigid  boundaries,  we  require  that  d  =  0  at  o  =  0, 1  which  simplifies 
the  thermodynamic  equation,  (3.2.14). 


In  Section  3.1,  we  considered  the  effects  of  curvature  and  beta  effect  on  the  planetary 
scale  flows  and  found  that  the  beta  effect  is  important  Now  we  consider  the  scale  for 
synoptic  flow  within  the  long  wave  environment.  We  anticipate  that  the  beta  effect  is  not  as 
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important  for  this  scale  of  flow.  Here,  we  examine  typical  synoptic  scales  for  the 
atmosphere  and  the  annulus  following  Chamey  (1948).  By  comparing  the  magnitudes  of 
various  terms  in  the  equations,  we  eliminate  terms  based  on  scaling  arguments.  Also,  we 
present  a  more  detailed  description  of  the  annulus  experiment. 

We  choose  the  following  scales  for  the  atmosphere  or  the  annulus: 

(1)  (x,y)*  [=]  L(x,y),  horizontal  length  scale, 

(2)  a*  [=]  D  o,  vertical  length  scale,  average  depth  of  the  fluid, 

(3)  (u,v)*  [=]  U(u,v),  horizontal  motion, 

(4)  w*  [=]  Uh/L(w),  vertical  motion,  from  scaling  the  continuity  equation, 

(5)  T'*  [=]  ATpj(T'),  scale  for  temperature  perturbations, 

(6)  Ts*  [=]  ATy(  Ts)  -  scale  for  the  imposed  temperature  gradient, 

where  *  represents  dimensional  quantities.  By  including  the  scale  in  these  equations,  all  of 
the  dependent  and  independent  variables  are  dimensionless.  If  we  combine  these  scales 
into  the  governing  equations  in  place  of  the  corresponding  variables,  we  can  non- 
dimensionalize  the  equations  by  dividing  through  by  the  scales  on  one  term.  Following 
this  procedure,  we  find  the  non-dimensional  governing  equations  are 


edu 

dt 


+  ^-  +  8  (a- 1)  Tv  =  0, 
dx  dx 


(3.3.1) 


+  U  +  3—  =  0, 

dt  dy  (3.3.2) 


tr*1-81*- 

V-  V  +  8~ — 2 — u^  =  0, 
do  (1  -  8tt)  dx 


(3.3.3) 


and 


(3.3.4) 
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where 


af+va-fi* 


dn 

a-u^a-l) 


dT, 


dz 


*«<>, 


<L  =-+v.y  +  5c^- 
dt  dt  ’ 


e  =  =  Rossby  number, 

LL* 


(3.3.5) 


V=M 
e  ’ 


S  =  =  Burger  number, 

and  tj  =  cos(2Kx)  which  represents  the  topography.  5  is  a  non-dimensional  mountain 
height.  The  Burger  number  may  be  thought  of  as  a  measure  of  the  stratification  or  a 
measure  of  the  Rossby  radius  of  deformation  to  the  horizontal  length  scale.  In  the 
remainder  of  this  work,  we  will  assume  that  is  not  a  function  of  z.  The  Rossby  radius 
of  deformation  gives  the  scale  for  which  we  expect  baroclinic  instability  to  occur. 

The  non-dimensional  parameter,  v,  is  a  measure  of  the  relative  strength  of  adiabatic 
cooling  versus  horizontal  temperature  advection.  When  we  apply  the  boundary  condition  o 
=  0  at  a  =  0,1,  (3.3.5)  will  be  applied  at  the  upper  and  lower  boundaries.  At  the  upper 
boundary,  the  last  term  on  the  left  hand  side  of  the  equation  is  zero  and  there  will  only  be 
temperature  advection.  At  the  lower  boundary,  v  becomes  important.  For  a  given 
stratification  and  Rossby  number,  the  magnitude  of  the  mountain  effect  is  dominated  by  the 
mountain  height  The  last  term  in  (3.3.5)  implies  that  any  zonal  flow,  u,  which  flows  over 
the  mountain  slope,  dq/dx,  will  rise  and  cool  adiabatically.  This  causes  an  imbalance  in 
(3.3.5).  To  overcome  the  imbalance,  a  southerly  flow  is  induced  which  advects  warm  air 
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from  the  south. 

The  mountain  in  this  problem  is  slowly  varying  in  the  x-direction  and  has  a  large 
horizontal  length  scale  compared  to  the  Rossby  radius  of  deformation  so  that  S  <  0(1). 

For  smaller  horizontal  length  scale  mountains,  S  ~  0(1).  Therefore,  in  this  model, 
baroclinic  instability  will  occur  for  scale  of  motions  where  L  is  approximately  Lp  and  S  ~ 
0(1).  This  argument  about  length  scales  suggests  that  the  mountain  will  force  planetary 
scale  flows  by  the  adiabatic  cooling  and  southerly  advection  of  temperature  but  should  have 
less  impact  on  the  growth  of  disturbances.  Since  baroclinic  eddies  propagate  within  the 
planetary  scale  flow,  the  large  scale  flow  will  act  as  a  wave  guide  for  the  smaller  scale 
eddies. 

The  magnitude  of  v  will  determine  how  large  the  mountain  appears  to  the  fluid.  If  v  is 
small,  due  to  any  combination  of  Burger  number,  Rossby  number  or  mountain  height,  the 
mountain  effect  is  small.  If  v  is  an  0(1)  quantity,  the  mountain  effect  will  become  very 
important. 

Typical  dimensional  and  non-dimensional  parameters  for  synoptic  scale  motions  in  the 


atmosphere 

;  are 

0) 

U=  10  ms'1 

(2) 

W  =  Uh/L  =  10'2 

(3) 

L  =  106  m 

(4) 

D  =  104  m 

(5) 

L/U  =  105  s 

(6) 

h  =  1(P  m 

(7) 

6  =  h/D  =  .1 

(8) 

e  =  .1 

(9) 

S  =  .25 
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(10)  Ld=  106  m 

(11)  v  =  .25 

(12)  N  =  2.5  x  10-3  s'1 

D  is  taken  as  a  typical  measure  of  the  depth  of  the  troposphere,  L  is  a  measure  of  the 
horizontal  length  of  the  phenomena.  It  is  important  to  distinguish  between  the  synoptic 
scale  and  the  planetary  scale  L.  U  is  an  average  value  of  the  wind  through  the  troposphere. 
W  can  be  estimated  from  the  thermodynamic  or  vorticity  equations,  but  in  this  problem  we 
find  the  scale  for  W  from  the  continuity  equation  and  it  is  related  to  the  mountain  height. 
When  the  fluid  interacts  with  the  mountain,  the  magnitude  of  W  will  be  larger  than 
predicted  by  the  thermodynamic  or  vorticity  equations.  We  have  taken  an  average  value  for 
the  mountain  height  to  be  1  km  yielding  8  =  0.1  which  is  the  same  order  as  the  Rossby 
number  for  synoptic  scale  motions.  As  we  will  see,  this  is  a  requirement  for  the  quasi- 
geostrophic  approximation  to  be  valid.  If  the  depth  of  the  mountain  were  much  larger  than 
the  Rossby  number,  the  leading  order  balance  would  not  be  geostrophic.  We  have 
assumed  the  stratification  is  a  linear  function  of  z  only  and  is  not  affected  by  the  mountain. 
A  value  of  S  =  0.25  means  that  the  length  scale  of  the  baroclinic  disturbances  is  about  one 
fourth  the  length  scale  of  observed  atmospheric  phenomena. 

Before  listing  parameter  values  for  the  annulus,  we  discuss  the  annulus  experiments  in 
more  detail.  The  interested  reader  should  refer  to  Pfeffer  et  al.  (1989)  or  Li  et  al.  (1986) 
for  a  complete  discussion  of  the  annulus,  experimental  procedures,  and  results. 

To  create  an  environment  that  simulates  the  Earth’s  atmosphere  in  a  controlled  way, 
Pfeffer  et  al.  (1989)  experiment  with  a  rotating  annulus.  Briefly,  the  annulus  is  a  rotating 
cylinder  that  is  filled  with  silicone  fluid.  It  is  heated  on  the  outer  wall  and  cooled  on  the 
inner  wall  so  that  a  radial  temperature  gradient  develops  and  a  baroclinic  flow  results.  The 
imposed  temperature  gradient  is  maintained  at  2°C.  A  two-wave  sinusoidal  topography 
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was  placed  on  the  bottom  of  the  cylinder.  For  the  experiments  of  interest  to  us,  the 
mountain  does  not  have  any  radial  structure  and  there  is  no  lid  in  contact  with  the  fluid. 
Figure  3.2  (a)  shows  the  top  view  of  the  annulus  dimensions.  The  radius  to  the  outer  wall, 
b  =  14.9  cm  and  to  the  inner  wall,  a  =  7.3  cm.  The  circumference  at  mid-radius  is  7t(a+b). 
An  azimuthal  view  of  the  annulus  is  shown  in  Figure  3.2  (b).  For  the  length  of  the 
channel,  we  take  the  circumference  at  mid-radius.  Therefore,  the  length  of  one  mountain  is 
half  the  channel  length.  We  choose  the  length  scale  for  the  problem  to  be  related  to  the 
mountain  by 

I  _  Lmountain 
2  K 

and  the  fundamental  zonal  wavenumber  is 

K  =  — 2eL —  = 

Lchannel 

and  the  fundamental  meridional  wavenumber  is 

.  _  7t  (a  +  b) 

“  4  (b  -  a)  ' 

Typical  dimensional  and  non-dimensional  parameters  for  the  annulus  are  given  as  follows: 

(1)  L  =  5.55  cm. 

(2)  D=  10  cm. 

(3)  h  =  1.75  cm. 

(4)  8  =  .175. 

Pfeffer  et  al.  (1989)  and  Li  et  al.  (1986)  ran  a  series  of  experiments  at  various  rotation 
rates.  The  remainder  of  the  parameters  of  the  problem  depend  on  the  rotation  rate.  Table 
3.1  shows  the  important  parameters  for  the  problem.  As  we  will  see,  for  the  0(8  *) 
approximation  to  be  valid,  the  ratio  of  e/8  must  be  approximately  0(1).  Table  3.1  shows 
this  quantity  is  less  than  0(1)  for  all  rotation  rates.  However,  we  will  show  in  Section  3.4 
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Figure  3.2  (a)  Top  view  of  the  annulus  where  b  is  the  outer  radius  and  a  is  the  inner 
radius,  and  (b)  side  view  of  the  annulus  dimensions  where  D  is  measured  from  the  surface 
z=0. 
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that  the  ratio  e/5  does  not  appear  in  the  quasi-geostrophic  potential  vorticity  equation. 
Rather,  v  =  S5/e  appears  in  the  thermodynamic  equation  at  the  lower  boundary.  This 
quantity  is  0(1)  for  all  rotation  rates.  Therefore,  we  expect  our  scaling  to  be  valid  for  most 
rotation  rates. 

Using  these  scale  parameters,  we  can  scale  the  governing  equations.  If  we  assume  that 
5  ~  O(e),  we  can  rewrite  the  problem  using  an  asymptotic  expansion  in  the  small 
parameter  5  such  that 

P  =  P0+  &  Pi  .... 

and 

r  =  To  +  8Tl  ... 

etc...  (3.3.6) 


Table  3.1  Typical  scale  parameters  for  the  annulus  as  a  function  of  rotation  rate. 
U  is  a  time  averaged  value  from  6.5  cm  at  mid-radius  are  listed.  L  =  5.55cm.  The 
topography  is  cos(2Kx). 


Q  (rad  s'1! 

0.628 

0.785 

1.05 

1.57 

2.09 

3.14 

U  (cm  s'1) 

0.32 

0.44 

0.29 

0.15 

0.11 

0.08 

Ld  (cm) 

3.62 

2.89 

2.16 

1.44 

1.09 

0.72 

e 

0.046 

0.050 

0.024 

0.009 

0.005 

0.002 

e/5 

0.262 

0.285 

0.137 

0.051 

0.029 

0.011 

S 

0.425 

0.271 

0.151 

0.067 

0.039 

0.017 

V 

1.62 

0.95 

1.10 

1.30 

1.37 

1.49 

Collecting  terms  of  the  same  order,  we  find  the  leading  order,  non-dimensional 
problem,  0(5°),  to  be 


-v 


o 


=  0, 


(3.3.7) 
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■*o 


+  ^°=0 

dy  ’ 


To, 


and 


where 


and 


(£)04+j<p°'> 


=  1. 


(3.3.8) 

(3.3.9) 

(3.3.10) 

(3.3.11) 


To  leading  order,  the  flow  is  geostrophic.  The  boundary  conditions  are  that  o  =  0  at  the 
upper  and  lower  boundaries. 

At  the  next  order,  0(8  ^),  we  find  the  correction  to  geostrophy  to  be 


and 


£ 

vsi 


dvp 

dt 


■M-l+^UO, 


V.V1+^.uo£=0, 


(3.3.12) 

(3.3.13) 

(3.3.14) 


where  the  ratio  e/8  is  to  be  viewed  as  an  0(1)  quantity  for  the  atmsophere.  However,  the 
magnitude  of  this  ratio  is  0(  10*  * )  for  the  annulus.  We  can  justify  keeping  this  term  in  spite 
of  the  scaling.  The  effect  of  topography  is  greatest  at  the  lowest  levels  of  the  fluid  but  is 
smallest  at  the  upper  level  where  advection  is  most  important.  In  (3.3.12)  when  a  =  0  at 
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the  bottom  of  the  fluid,  the  balance  may  be  between  the  last  three  terms  on  the  left  hand  side 
of  this  equation.  However,  when  a  >  0,  the  last  term  on  the  left  hand  side  of  (3.3.12)  is 
small  and  now  the  balance  must  be  between  the  first  three  terms.  For  this  order,  we 
included  only  those  equations  required  for  the  following  derivations. 


Next,  we  can  form  the  quasi-geostrophic  potential  voracity  equation  using  the  0(6°) 
and  0(6  * )  equations  outlined  in  Section  3.3.  The  static  temperature  was  taken  as  a 
function  of  the  z-coordinate  only  when  we  transformed  the  equations  to  the  sigma 
coordinate  system.  Physically,  this  makes  good  sense  since  in  the  absence  of  motion,  the 
stratification  is  a  function  of  z  and  does  not  vary  along  the  new  vertical  coordinate.  This 
assumption  simplifies  the  final  form  of  the  equations.  If  we  take  the  static  temperature  to 
be  a  linear  function  of  the  vertical  coordinate,  the  Brunt  Vaisala  frequency  becomes  a 
constant.  We  can  form  the  vertical  voracity  equation  by  taking  5(3.3. 12)/5y  minus 
5(3.3. 13)/5x.  Next,  we  can  solve  for  a  in  (3.3.1 1)  and  substitute  for  vQ  using  (3.3.7),  u0 
using  (3.3.8),  and  TQ  using  (3.3.9)  and  obtain  the  expression 


i|  +  J(Po-> 


5o\  5y  /  5y\  do  I 
(3.4.1) 


where  the  last  term  on  the  right  hand  side  is  the  thermal  wind  balance  and  cancels.  The 
second  term  on  the  right  hand  side  is  also  zero  but  this  is  true  only  since  we  took  T$  to  be  a 
function  of  z.  There  is  a  vertical  derivative  of  the  static  temperature  implied  in  the  second 
term  but  we  took  this  term  to  be  one.  If  we  combine  the  vertical  vorticity  equation  and 
(3.4.1),  the  quasi-geostrophic  potential  vorticity  equation  can  be  written  as 
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v2p°+srfra 

s  °<*ll  (3.4.2) 

In  the  final  form  of  the  vorticity  equation,  only  the  Burger  number  appears  as  a  non- 
dimensional  parameter  and  the  ratio  of  8/e  canceled  out  of  the  problem. 

The  thermodynamic  equation  with  the  boundary  condition  a  =  0  at  the  upper  and  lower 
surfaces  may  be  written  as 


^  +  J(PC.  ) 


d_ 

dt 


+  (a- 1)  v 


apQ 

dy  dx 


=  0. 


(3.4.3) 


The  last  term  on  the  left  hand  side  of  the  equation  only  has  a  value  when  a  =  0,  at  the  lower 
boundary.  It  represents  the  interaction  of  the  low  level  zonal  flow  with  the  mountain.  The 
non-dimensional  number  v  =  (S8/e)  controls  the  magnitude  of  this  term  which  physically 
describes  adiabatic  heating  and  cooling  due  to  the  flow  over  the  topography.  If  the  non- 
dimensional  mountain  height,  5,  is  small  compared  to  the  Rossby  number,  e,  then  the 
mountain  will  not  play  a  major  role  in  the  problem.  If  v  ~  0(1),  then  this  term  is  important 
in  balancing  the  full  thermodynamic  equation.  If  the  mountain  height  is  very  large 
compared  to  the  Rossby  number  so  that  v  >  0(1),  then  the  approximations  used  to  derive 
the  quasi-geostrophic  equations  are  invalid.  Malguzzi  (1987)  showed  that  the  quasi- 
geostrophic  approximation  is  invalid  for  steep  mountains.  The  geostrophic  and  hydrostatic 
approximations  allow  us  to  write  the  pressure,  P0,  in  terms  of  the  streamfunction,  \jr  =  P0, 
so  that  (3.4.2)  and  (3.4.3)  become 


(3.4.4) 


a_ 

dt 


(3.4.5) 
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Equations  (3.4.4)  and  (3.4.5)  constitute  the  governing  equations  for  the  remainder  of  this 
work. 


_ Linearized  Quasi-Geostrophic  Model 


We  are  primarily  interested  in  testing  the  stability  of  basic  states  representative  of  flow 
over  topography.  To  that  end,  we  seek  to  linearize  the  equations  about  a  basic  state  and  so 
we  split  the  total  streamfunction  into  basic  state  and  perturbation  components  as  follows 


\|/(x,y,<T,t)  =  ^(x.y.o)  +  (p(x,y,o,t)  (3.5.1) 

Here,  'P  represents  the  steady  basic  state  which  in  this  problem  is  a  function  of  the  x- 
coordinate  as  well  as  the  y-  and  a-coordinates.  Substituting  for  \j/  in  (3.4.4)  with  (3.5. 1 ), 
we  find 

+  J(^,n)  +  J(^,q)  +  J((p,n)  +  J(cp,q)  =  0, 
dt  (3.5.2) 

(a)  (b)  (c)  (d)  (e) 

where  the  perturbation  potential  vorticity,  q,  is  defined  as 


q  =  V2<p  + 


1  d2<P 
S  do2 


and  the  basic  state  potential  vorticity,  n,  is  defined  as 

s  0c2 

In  (3.5.2),  term  (a)  represents  the  local  rate  of  change  of  the  perturbation  potential  vorticity. 
Term  (b)  represents  the  basic  state  advection  of  the  basic  state  potential  vorticity.  Term  (c) 
represents  represents  the  basic  state  advection  of  the  perturbation  potential  vorticity.  Term 
(d)  represents  the  perturbation  advection  of  the  basic  state  potential  vorticity.  Term  (e) 
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represents  the  nonlinear  perturbation  advection  of  the  perturbation  potential  vorticity. 

Substituting  for  iy  in  (3.4.5)  with  (3.5.1),  the  thermodynamic  boundary  conditions 
may  be  written  as 


a/a<p\ 

at\ao) 

+  jj'p  ^.1 
T’da/ 

+jh!) 

•tel 

+ 

=  0 

(3.5.3) 

(a) 

(b) 

(d) 

(e) 

(0 

at  0  =  1 ,  and 

33*' 

U,  d^1 

I’doj 

+  jKvti) 

♦te 

) +  V’te 

+  VTlj  +  j| 

1*55) 

=  0 

(a) 

(b) 

(c) 

(d) 

(e) 

<0 

(3.5.4) 


at  a  =  1.  Equations  (3.5.3)  and  (3.5.4)  are  the  thermodynamic  boundary  conditions  at  the 
upper  and  lower  boundaries.  Equation  (3.5.4)  contains  all  of  the  terms  that  appear  in 
(3.5.3)  and  the  additional  mountain  terms  for  the  lower  boundary.  Term  (a)  represents  the 
local  rate  of  change  of  the  temperature.  Term  (b)  represents  basic  state  advection  of  the 
basic  state  temperature.  Term  (c)  represents  advection  along  the  mountain  slope.  Term  (d) 
represents  basic  state  advection  of  perturbation  temperature.  Term  (e)  represents 
perturbation  advection  of  the  basic  state  temperature  and  the  flow  along  the  mountain  slope. 
Term  (f)  represents  the  nonlinear  perturbation  advection  of  the  perturbation  temperature. 

For  the  remainder  of  this  problem,  we  select  the  linearized  set  of  equations  with 
respect  to  the  perturbation  and  neglect  the  last  term  in  (3.5.2),  (3.5.3),  and  (3.5.4).  In 
Section  3.6,  we  will  find  an  x-dependent  basic  state  which  exactly  satisfies  all  of  the  terms 
in  each  of  these  equations  that  do  not  contain  the  perturbation.  In  (3.5.2)  and  (3.5.3),  term 
(b)  will  be  exactly  satisfied  and  equal  to  zero  by  our  choice  of  basic  state.  In  (  3.5.4), 
terms  (b)  and  (c)  will  sum  to  zero  by  our  choice  of  basic  state. 
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3.6  Wavv  Basic  State 

In  Chapter  2,  we  discussed  the  climatological  mean  and  low  frequency  -m  lability  in  the 
Northern  Hemisphere  winter  circulation  and  the  fact  that  low  frequency  variability  accounts 
for  the  majority  of  the  climatological  distribution.  If  we  assume  these  waves  are  standing, 
we  can  fin  teady  state  solutions  to  the  quasi-geostrophic  potential  vorticity  and 
thermodynamic  equations  at  the  upper  and  lower  boundaries 

The  basic  state  we  choose  must  be  an  exact  solution  to  the  quasi-geostrophic  potential 
vorticity  equation  as  well  as  the  thermodynamic  boundary  conditions.  One  choice  is  to 
force  the  flow  at  the  lower  boundary  to  go  to  zero.  This  is  a  poor  choice  since  there  would 
be  no  interaction  of  the  basic  state  flow  with  the  mountain  slope.  Then,  the  interaction  with 
the  mountain  must  come  into  the  problem  at  the  order  of  the  perturbation.  A  better  choice  is 
to  allow  for  flow  at  the  bottom  boundary  and  let  it  interact  with  the  mountain  but  this 
forces  us  to  choose  an  x-dependent  basic  state. 

Held  and  Ting  (1990)  stress  the  importance  of  the  low  level  flow  in  forcing  stationary 
waves.  Using  a  numerical  model  with  topography,  they  show  that  the  amplitude  of  the 
forced  stationary  wave  grows  as  the  flow  at  the  surface  increases.  In  the  atmosphere,  the 
constant  flow  at  the  surface  results  from  momentum  flux  convergence  into  the  mid-latitudes 
from  the  tropics  and  in  general  the  flow  is  westerly.  This  analytical  model  does  not  have 
side  wall  boundary  conditions  that  can  transport  momentum  from  the  tropics.  Therefore, 
we  must  parameterize  the  constant  zonal  flow. 

There  have  been  many  studies  of  the  stability  properties  of  zonal  basic  flows  but  there 
are  comparatively  fewer  stability  studies  using  wavy  basic  states.  Grotjahn  (1984)  wrote  a 
review  article  on  wavy  basic  states  used  in  linear,  theoretical  models  of  baroclinic 
instability.  The  majority  of  the  works  used  two-layer,  beta  plane  models.  Niehaus  (1980) 
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used  a  wavy  basic  state,  and  a  continuous  f-plane  to  do  a  linear  instability  analysis  of  a 
baroclinic  flow.  The  choice  of  basic  state  was  arbitrary.  It  was  not  forced  by  topography 
or  heating  but  rather,  it  was  chosen  to  match  the  observed  conditions  in  the  Northern 
Hemisphere  winter.  She  compared  wavy  basic  states  and  skewed,  wavy  basic  states  and 
found  that  the  tilt  in  the  basic  state  allows  for  a  natural  interaction  between  the  forced  wave 
and  the  instabilities.  Without  the  skewness  in  the  basic  state,  the  fastest-growing 
instabilities  had  zero  meridional  wavenumber.  For  skewed  basic  states,  the  fastest  growing 
instabilities  were  coupled  to  and  modulated  by  the  standing  wave.  Buzzi  et  al.  (1984) 
studied  the  stability  of  a  wavy  basic  state  in  a  two-layer  beta-plane  model  and  a  continuous 
f-plane  model.  In  their  linear  stability  analysis,  the  instability  had  the  same  wavenumber  as 
the  topography.  They  found  topographic  instability  occurs  for  long  waves;  long  waves 
which  are  beyond  the  cutoff  for  normal,  traveling  baroclinic  waves.  They  compare  the 
growth  rates  for  wavy  and  non-wavy  basic  states  in  both  models  and  found  growth  rates 
are  slightly  decreased  for  the  case  of  a  wavy  basic  state. 

Derome  (1984)  shows  one  procedure  to  find  finite  amplitude ,  steady-state,  x-dependent 
basic  states  which  represent  exact  solutions  to  the  nonlinear  vorticity  equation  and 
boundary  conditions.  The  importance  of  his  work  is  that  under  certain  circumstances, 
rather  complex  solutions  to  the  linearized  potential  vorticity  equation  are  also  finite- 
amplitude  solutions  to  the  equations.  His  procedure  is  limited  to  inviscid  problems. 

In  this  section,  we  follow  Derome ’s  procedure  to  find  an  x-dependent,  steady-state, 
baroclinic  basic  state.  We  assume  the  x-dependent  basic  state  can  be  separated  into  one  part 
that  has  linear  vertical  shear  and  an  x-dependent  part  such  that 

'F(x,y,o)  ='FB(y,o)  +  <|fe(x,o) 

’?B(y.a)  =  -y(u0  +  Mo+^)). 


where 


(3.6.1) 
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We  used  the  subscript  B  to  denote  basic  state  quantities.  Eady  (1949)  used  a  basic  state 
with  linear  vertical  shear  to  show  the  simplest  form  of  baroclinic  instability.  Here  X 
represents  the  vertical  shear  and  U0  is  a  constant  zonal  flow.  For  convenience,  we 
redefined  the  vertical  coordinate  so  that  the  lower  boundary  is  at  o  =  -0.5  and  the  upper 
boundary  is  at  o  =  +0.5.  If  we  substitute  this  basic  state  into  the  quasi -geostrophic 
potential  vorticity  equation  and  boundary  conditions,  we  can  rewrite  the  equations  as 


1  anB 


u  dy 


<t>B  + 


=  o, 


(3.6.2) 

(3.6.3) 


where 


_ fek  +  vn)  + id>n.^-  +  vn) 


U  =  - 


(3.6.4) 


and 


qB  =  v  + 


l  d2<te 
S  ’ 


nB  =  vVB  + 


i  a2^B 

s  ao2  ' 


Equation  (3.6.2)  is  the  potential  vorticity  equation  and  (3.6.3)  and  (3.6.4)  are  the  upper 
and  lower  boundary  conditions.  We  note  that  if  the  first  term  on  the  left  hand  side  of 
(3.6.2)  -  (3.6.4)  is  zero,  then  the  second  term  also  must  vanish.  This  implies  that 


qB  =  - 


i  anB 
u  dy 


(3.6.5) 


is  one  of  many  exact  solutions  of  the  basic  state  potential  vorticity  equation.  To  simplify 
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further,  we  chose  the  Eady  basic  state  so  that  rig  =  0.  It  follows  from  (3.6.5)  that  qg  =  0. 
We  could  have  taken  qg  to  be  some  other  function  of  the  vertical  coordinate.  A  better 
choice  would  have  included  the  (3  effect  A  general  solution  for  the  streamfunction  is 


where 


4>b(x,0)  =  vFg(a)  cos(2Kx) 

Fg(o)  =  Ag  cosh(y(c  -  ^))  +  bb  sinh(y(a  -  ^))  f 
Y  =  2KVS\ 


x-a  cross  section  of  the  x-dependent  part  of  the  basic  state.  The  mountain  ridge  is  at  0  and 
x  rad.  The  high  pressure  is  over  the  mountain  ridge.  Since  the  flow  is  inviscid  and  steady, 
there  is  no  tilt  of  the  basic  state  with  height.  If  the  flow  were  viscous,  there  would  have 
been  a  tilt  with  height  because  the  steady  flow  would  extract  energy  from  the  basic  state  to 


Figure  3.3  Basic  state  streamfunction  where  (a)  is  an  x-y  cross  section  at  the  top,  (b)  is 
and  x-y  cross  section  at  the  bottom,  (c)  is  an  x-c  cross  at  mid-radius,  and  (d)  is  an  x-ct 


cross  section  of  the  wavy  part  of  the  basic  state  at  mid-radius.  Mountain  ridges  are  at  0  and 
it  rad.  S  =  0.25,  U0=  0.1,  5  =  0.1,  and  X  =  1. 
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overcome  dissipation.  These  figures  may  be  compared  to  the  January  climatology  in 
Figures  2.2,  2.5,  and  2.8. 

Figure  3.4  shows  the  amplitudes  of  the  coefficients  Ag  and  Bg  in  (3.6.7)  as  a  function 
of  the  basic  flow,  UQ.  For  negative  values,  the  basic  state  is  very  sensitive  to  the 
magnitude  of  UD  especially  as  U0  approaches -0.5.  This  happens  because  the 
denominator  of  A  which  is  proportional  to  (U0  +  hi)  becomes  very  small.  Thus,  as  UQ 
approaches  -0.5,  the  denominator  gets  very  small  and  the  amplitude  of  this  term  becomes 
very  large.  Negative  values  of  UQ  have  a  significant  impact  on  the  growth  rates  and 
frequencies  of  instabilities  on  this  basic  state  as  we  will  show  in  Chapter  4. 


tTnTiTTTTn  ITral  iTm  pi  ol  n  sWTtTi  ■  DWiTTB  OTT- 


Unlike  typical  baroclinic  and  barotropic  instability  where  the  unstable  waves  propagate 
and  grow,  topographic  instability  is  characterized  by  waves  which  grow  but  do  no 
propagate.  Several  papers  have  been  written  on  topographic  instability.  Chamey  and 
DeVore  (1979,  hereafter  CD)  showed  topographic  instability  occurs  in  a  barotropic,  beta- 
plane  channel  model  with  dissipation  and  topography.  Chamey  and  Straus  (1980,  hereafter 
CS)  showed  topographic  instability  occurs  in  a  baroclinic,  two-layer,  beta-plane  channel 
model  with  dissipation  and  topography.  Buzzi  etal.  (1984)  showed  topographic  instability 
occurs  in  a  baroclinic,  two-layer,  beta-plane  channel  model  with  dissipation,  zonally 
symmetric  and  asymmetric  basic  states,  and  topography  and  also  in  a  continuous,  f-plane 
channel  model  with  dissipation,  zonally  symmetric  and  asymmetric  basic  states,  and 
topography.  Yoden  and  Mukougawa  (1983)  used  a  baroclinic  two-layer  and  20-layer, 
low-order  model  to  investigate  flow  over  topography. 

Kallen  (1983)  found  topographic  instability  in  a  baroclinic  model  that  included  a  direct 


47 


momentum  forcing  in  the  vorticity  equations  in  the  lower  level  of  his  two  layer  model.  He 
showed  that  topographic  instability  does  not  occur  in  baroclinic  models  unless  a  momentum 
source  is  included  and  pointed  to  some  inconsistencies  in  the  work  of  CS.  In  our 
baroclinic  model,  the  momentum  source  is  the  constant  zonal  flow,  UQ. 

Revell  and  Hoskins  (1984,  hereafter  RH)  critically  reviewed  CD’s  barotropic  model  and 
CS’s  baroclinic  model.  In  their  work,  RH  wanted  to  verify  if  topographic  instability  was 
different  from  baroclinic  and  barotropic  instability.  They  present  a  thorough  review  of 
topographic  instability  in  a  barotropic  and  a  baroclinic  model.  They  found  that  topographic 
instability  does  exist  in  a  barotropic  model.  However,  in  baroclinic  models,  the 
topographic  instability  appeared  in  simple,  low  order  models  but  did  not  appear  in  more 
complete  baroclinic  models  where  classical  baroclinic  instability  dominated  at  all 
wavenumbers.  As  RH  pointed  out ,  it  would  be  difficult  to  verify  topographic  instability 
using  real  atmospheric  data,  however  it  does  appear  that  the  instabilities  of  the  linearized 
problem  play  an  important  role  in  highly  truncated,  simple  models.  In  this  work,  we  will 
place  more  emphasis  on  trying  to  understand  how  baroclinic  instability  is  modified  rather 
than  on  topographic  instability. 

Buzzi  et  al.  (1984)  outline  the  process  that  leads  to  topographic  instability  in  their 
models.  The  physical  explanation  for  the  instability  starts  with  the  form-drag  mechanism. 
The  topography  exerts  drag  on  the  fluid  due  to  surface  pressure  components  out  of  phase 
with  the  topography.  Drag  leads  to  momentum  exchanges  between  the  wavy  disturbances 
and  the  mean  zonal  flow.  This  implies  there  is  a  form-drag  coupling  between  wavy 
disturbances  and  the  mean  zonal  flow  caused  by  topography.  Form-drag  coupling 
accelerates  the  mean  zonal  flow.  We  can  derive  the  form-drag  equation  from  the  x- 
component  of  the  momentum  equation.  Restating  the  0(8  * )  x-component  of  the 
momentum  equation,  where  sigma  now  varies  from  -0.5  to  0.5 
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<°-i>!To=0' 


(3.7.1) 


We  can  rewrite  this  equation  in  flux  form  as 

(*)(£  +  *  <“oV-o) -(uo  +  vo)V.V-oj.v1+^i  +  (a.i)|To  =  0. 

(3.7.2) 

If  we  substitute  the  leading  order  continuity  equation,  (3.3.9),  and  hydrostatic  equation, 
(3.3.8),  and  integrate  over  volume,  (3.7.2)  may  be  written  as 

L  lf)(^  +  *  ^  • V1  +  ^  +  <0‘  2^T°! dvo1  =  °- 

(a)  (b)  (c)  (3.7.3) 


In  our  model,  we  chose  the  streamfunction  to  be  periodic  in  x  and  y.  Therefore,  terms 
(a)  and  (c)  are  zero.  Term  (b)  is  zero  since  there  is  no  net  north-south  mass  flux.  Equation 

(3.7.3)  reduces  to  the  form 


1  Mdy 

2  ox  Bo 


dvol  =  0. 


(3.7.4) 


In  (3.7.4),  the  first  term  on  the  left  hand  side  represents  the  acceleration  of  the  zonal  flow 
due  and  the  second  term  represents  a  mountain  torque.  In  Chapter  4,  we  will  substitute  for 
\\f  using  our  assumed  solutions  and  complete  the  integration  of  (3.7.4).  We  will  use  this 
equation  to  investigate  topographic  instability  and  show  that  topographic  instability  exists 
for  perturbations  with  the  same  wavenumber  as  the  topography  and  no  y-structure.  Based 
on  the  work  of  Revell  and  Hoskins  (1984)  and  Kallen  (1983),  we  do  not  expect 
topographic  instability  unless  there  is  flow  at  the  surface.  In  Chapter  4  we  solve  the 
governing  equations  and  include  a  constant  zonal  flow  U0. 


CHAPTER  4 


ANALYTICAL  AND  NUMERICAL  SOLUTIONS 

In  Chapter  3,  we  developed  a  mathematical  model  for  large  scale  and  synoptic  scale, 
baroclinic  flow  over  two-wave  bottom  topography.  We  chose  a  very  simple  model  to  focus 
on  two  important  features  of  the  flow;  the  large  scale  flow  interacting  with  the  mountain  to 
force  a  wavy  basic  state  and  the  synoptic  scale  flow  interacting  with  the  mountain  and  the 
basic  state  to  form  storm  track  regions  and  preferred  regions  for  growth.  We  neglected  the 
beta  effect,  friction,  and  curvature  that  would  make  our  model  more  correct  but  also  more 
complicated. 

Many  researchers  have  completed  linear  stability  analyses  and  studied  how  orography 
modifies  baroclinic  instability.  Kallen  (1983),  Yoden  and  Mukougawa  (1983),  Speranza 
etal.  (1985),  Buzzi  and  Speranza  (1986),  and  Buzzi  etal.  (1987),  Reinhold  (1990)  have 
studied  the  problem  using  two-layer,  baroclinic,  beta-plane  models  with  topography  and 
dissipation.  Yoden  and  Mukougawa  (1983)  studied  baroclinic  instability  using  a  20-layer 
sigma  coordinate,  beta-plane  model  with  topography  and  dissipation.  Niehaus  (1980)  and 
Speranza  etal.  (1985)  studied  modified  baroclinic  instability  using  a  continuous,  f-plane 
model  with  topography  and  dissipation.  Malguzzi  et  al.  (1987)  studied  the  effect  of  finite 
height  topography  on  non-geostrophic  baroclinic  instability  and  discussed  the  implications 
for  simple  theories  of  baroclinic  instability.  DeSzoeke  (1983)  studied  the  destabilizing 
effect  topography  has  on  otherwise  stable  waves.  Buzzi  et  al.  (1987)  examined  the 
interaction  of  a  baroclinic  eddy  with  topography  using  a  linear,  two-layer  model.  They  use 
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only  zonal  basic  flows.  Therefore,  they  are  studying  how  the  perturbation  interacts  with 
topography.  Buzzi  etal.(  1984)  studied  the  effects  that  topography  has  on  the  growth  rates 
of  unstable  modes.  They  find  that  the  growth  rates  are  decreased,  especially  at  higher 
wavenumbers.  Most  of  the  other  theoretical  works  focus  on  this  same  aspect. 

In  this  chapter,  we  do  a  linear  stability  analysis  using  an  analytical  and  a  numerical 
technique.  The  analytical  model  is  an  extension  of  the  Eady  (1949)  model  with  a  correction 
due  to  topography.  The  solution  procedure  involves  an  asymptotic  expansion  in  the  small 
parameter  v  =  S8/e  in  the  limit  as  v  ->  0,  and  S  and  e  are  fixed.  The  numerical  solution 
also  incorporates  the  Eady  model  but  we  solve  the  problem  by  transforming  the  full  set  of 
equations  to  spectral  space  and  then  perform  an  eigenanalysis.  Finally,  we  consider  the 
topographic  instability  problem  where  the  perturbation  grows  but  does  not  propagate.  We 
begin  by  reviewing  the  Eady  model. 


4.1  Eady  Model 

One  of  the  simplest  models  of  baroclinic  instability  was  introduced  by  Eady  (1949).  In 
this  model  the  north-south  temperature  gradient  is  simulated  by  the  vertical  shear  of  the 
horizontal  wind  through  the  thermal  wind  relationship.  The  vertical  shear  is  taken  to  be 
linear.  The  flow  is  confined  between  two  rigid  boundaries  at  the  top  and  bottom.  The 
stratification  and  density  are  taken  as  constant  The  streamfunction  may  be  written  as 

\|r  =  -yaX+<p(x,y,a,t)  (4.1.1) 

where  the  first  term  on  the  right  hand  side  represents  the  basic  state  and  the  second 
represents  the  perturbation.  For  this  choice  of  basic  state,  n  =  0  and  (3.5.2)  reduces  to 
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with  boundary  conditions 

[da/  (4.1.3) 


<L  +  uT.B— 
at  ax 


where 

UT  =  -  ^  =  U0  +  X, 

ay 

U  "  9y  -  U0* 

and  the  superscript  T  is  for  the  top  boundary  and  B  is  for  the  bottom  boundary.  Solutions 
for  (4.1.2)  are  of  the  form 

<p(x,y,c,t)  =  F0(o)  ei(Knx  +  /my  -  ax)  +  CC,  (4  Y  4) 


where  CC  denotes  complex  conjugate.  Here,  F0(o)  represents  an  undetermined  function 
of  the  vertical  coordinate,  K  is  the  fundamental  zonal  wavenumber,  n  is  the  zonal 
wavenumber,  /  is  the  fundamental  meridional  wavenumber,  and  m  is  the  meridional 
wavenumber.  We  take  m  =  1  for  the  remainder  of  this  section.  Substituting  (4.1.4)  into 
(4.1.2),  we  find 


where 


(uT,BKn  -  co) 


92Fr 


do2 


-  4|1$  F0 


=  0 


4nJ  =  (K2n2  +  /2) 


(4.1.5) 


The  boundary  condition  at  a  =  +0.5  is 

(-co  +  KnUT)^-\Kn  =  0, 
do 

and  the  boundary  condition  at  o  =  -0.5 

(-O)  +  KnUB)^--XKn  =  0. 
do 


(4.1.6) 


(4.1.7) 


A  general  set  of  discrete  solutions  to  (4.1.5)  is  given  by 
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Fo(0)  =  Ancosh(  2p,n(o  -  i-))  +  Bnsinh(  2\ln(a  - 1))  (4  l  g) 

Applying  the  boundary  conditions,  we  obtain  the  following  set  of  algebraic  equations  for 
AnandBn 


-XKn 

(-to  +  UTKn)2Hn 

An 

(XKn  +  (-to  +  ITBKn)2|inTanh(2(ln)) 

((-to  +  UBKn)2p.n  +  X.KnTanh(2(in)) 

.  ®n  . 

=  0. 


(4.1.9) 

We  can  find  nontrivial  solutions  for  An  and  Bn  if  the  determinant  for  the  coefficients 
vanishes.  Solving  the  determinant,  we  find  the  dispersion  relation  for  the  Eady  model  may 
be  written  as 


_CQ_  = 

Kn 


=  U0  +  \  ±  a[1  -  -Mcoth(2lin)  -  -M 
°  2  14  S  2|Xn /] 


,i_L 

2 


(4.1.10) 

Figure  4.1  (a)  shows  a  plot  of  growth  rate  versus  wavenumber  for  S  ranging  from  0.1  - 
0.4.  For  S  =  0.25,  the  fastest  growing  mode  is  wavenumber  six.  Figure  4.1  (b)  shows  a 
plot  of  frequency  versus  wavenumber  for  S  ranging  from  0.1  -  0.4.  As  the  stratification 
increases,  the  growth  rate  decreases  and  the  frequency  increases.  For  short  waves,  there 
are  two  neutral  modes  which  have  different  frequencies.  For  long  waves,  beyond  the 
shortwave  cutoff,  there  are  two  roots  which  are  complex  conjugates.  One  of  the  roots  is 
growing  and  the  other  is  decaying.  Only  the  growing  root  is  shown  in  Figure  4.1  (a). 
These  waves  propagate  with  phase  speed  (UD  +  X/2). 

The  Burger  number,  S,  relates  the  Rossby  radius  of  deformation,  Lj),  to  the  length 
scale  of  the  phenomena  by 

\2 

(4.1.11) 


As  discussed  by  Pedlosky  (1987),  for  S  =  0.25,  the  dimensional  wavelength  is  given  by 
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Figure  4.1  (a)  Growth  rate  as  a  function  of  wavenumber  for  the  Eady  model.  S  ranges 
from  0.1  -0.4.  S  =  0.25  is  labeled,  (b)  Frequency  for  the  positive  root  in  (4.1.10)  and  (c) 
for  the  negative  root  in  (4.1.10)  as  a  function  of  wavenumber  for  the  same  range  of  S. 
Arrows  indicate  the  direction  of  increasing  stratification. 
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A,-27tL_4”LD 

Kn  Kn  '  (4.1.12) 

If  we  take  Lq  to  be  lO^km  for  the  atmosphere,  then  (4.1.12)  gives  a  wavelength  of  the 
order  of  4000  km  that  will  dominate  the  fluctuations  in  the  atmosphere.  This  is  in 
agreement  with  the  observed  synoptic  scale  disturbances. 

Certain  features  of  the  Eady  model  are  unrealistic  and  cannot  be  applied  to  the 
atmosphere  although  the  model  has  been  applied  successfully  to  the  annulus.  The  most 
important  feature  missing  from  the  model  is  the  planetary  vorticity  gradient.  Chamey  and 
Stem  (1962)  proved  there  can  be  no  growing  disturbances  in  a  basic  flow  unless  the  north- 
south  gradient  of  basic  state  potential  vorticity  changes  sign  somewhere  in  the  fluid.  In  the 
Eady  model,  the  potential  vorticity  is  constant  throughout  the  fluid.  Therefore,  there  is  no 
north-south  basic  state  potential  vorticity  gradient.  However,  Bretherton  (1966)  showed 
that  the  thermodynamic  boundary  conditions  are  equivalent  to  sheets  of  potential  vorticity 
just  above  and  below  the  fluid  boundaries.  Thus  there  are  two  infinitesimally  thin  sheets  of 
potential  vorticity  within  the  fluid  and  so  a  basic  state  potential  vorticity  gradient  exists. 

The  gradients  are  equal  in  magnitude  but  opposite  in  sign  at  the  top  and  bottom.  Therefore, 
we  need  an  upper  boundary  in  the  Eady  model  for  baroclinic  instability  to  develop. 

4.2  Analytical  Solutions  and  Results 

In  Section  4.1  we  reviewed  the  Eady  model.  In  this  section  we  will  use  the  Eady 
model  to  obtain  an  analytical  solution  for  the  problem  with  topography.  We  begin  by 
substituting  the  wavy  basic  state  we  derived  in  Section  3.6  into  (3.5.2)  -  (3.5.4).  We  find 
all  of  the  terms  involving  the  mountain  are  scaled  by  v.  Next,  we  assume  v  is  small  so  that 
we  can  rewrite  the  equations  using  an  asymptotic  expansion  in  the  small  parameter  v  = 
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(S8/e).  To  leading  order,  0(v°),  we  find  the  problem  is  the  Eady  model.  At  the  next 
order,  0(v  *),  we  bring  in  the  correction  due  to  the  mountain.  We  chose  the  model  to  be 
periodic  in  x  and  y.  There  is  a  two-wave  bottom  topography  and  a  rigid  upper  lid. 


Substituting  the  wavy  basic  state  into  (3.5.2)  gives 
— -  +  JOF  ,q)  =  0. 

From  (3.5.3),  the  boundary  condition  at  a  =  +0.5  becomes 


0404S)-* 


(4.2.1) 


(4.2.2) 


and  from  (3.5.4),  the  boundary  condition  at  cr  =  -0.5  becomes 


0+045 


+  VT|  =  0. 


(4.2.3) 


If  for  simplicity,  we  take  q  =  0  at  time  t  =  0,  then  q  =  0  for  all  time  and  (4.2. 1 )  becomes 


(4.2.4) 


Substituting  (3.6.10)  into  (4.2.2)  and  (4.2.3),  we  can  rewrite  the  boundary  conditions  as 

!  +  UT’BM!)  -^-v^KsinaRx))^  -^  +  (O-l)]^0. 

(a)  (b)  (c)  (4.2.5) 

where  =  (UD  +  X)  and  U®  =  UQ.  In  (4.2.5),  (a)  is  the  basic  state  meridional  advection 
of  the  meridional  perturbation  temperature  gradient,  (b)  is  the  zonal  perturbation  advection 
of  the  zonal  basic  state  temperature  gradient,  and  (c)  is  the  zonal  perturbation  flowing  over 
the  mountain  slope.  Recall  that  the  expression  for  Fg  is  given  by 


FB(o)  =  Ab  cosh  (y(a  -  ^ ) )  +  BBsinh  (y(c  -  ^)) 


(3.6.7) 
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and  Ag  and  Bg  are  given  by  (3.6.8)  and  (3.6.9).  If  we  assume  solutions  of  the  form 

<p(x,y,o,t)  =  x(x,a)ei(/y  -  <**)  +  CC,  (4.2.6) 

(4.2.4)  may  be  written  as 


dx2  ^dd2 


and  (4.2.5)  may  be  written  as 


-io)+UT3^- 

dx 


dy  dy 

~  v(i/2Ksin(2Kx)) 

do  dx 


Fd—  -  ^F-B.  +  (o .  1) 
Bdo  do  1  2 


(4.2.7) 

X  =  0. 

(4.2.8) 


4.2.2  Asymptotic  Expansions  In  Orders  of  v 

We  assume  that  we  can  rewrite  the  variables  of  (4.2.7)  and  (4.2.8)  in  terms  of  the 
small  parameter  v  =  (S6/e)  such  that 


X  =  Xo  +  vlXl  +•••. 


and 


o)  =  (ob  +  v‘co1  +  ....  (4.2.9) 

We  can  substitute  (4.2.9)  into  (4.2.7)  and  (4.2.8)  and  collect  terms  by  order.  The  0(v°) 


problem  becomes 


Idx2 


IV  I  la2Xo_0 

1X0  s  302  -° 


(4.2.10) 


and 


(4.2.11) 


We  find  the  0(v  * )  problem  to  be 


(iMxl+©=° 


and 


(4.2.12) 
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'  UT'B|r]  (If)  -  *ff '  -  <^Sin(2Kx))[f 


The  leading  order  problem  is  the  Eady  problem  that  was  solved  in  Section  4.1.  For  this 
problem,  we  assume  solutions  of  the  form 

Xo  =  F„(o)eiK"*  +  CC  (4.2.14) 

and  substitute  into  (4.2.10)  and  (4.2.1 1)  to  obtain  the  Eady  dispersion  relation  given  by 
(4.1.10).  The  growth  rates  and  frequencies  are  shown  in  Figure  4.1. 

The  O(v^)  equations  also  involve  leading  order  quantities.  Therefore,  we  substitute 
the  leading  order  solution,  (4.2.14)  into  (4.2.13)  and  obtain 

-in*,  +  UT>Bl-  j^-J .  =  iojj^eiKnx  +  (/K)  (eiK(n+2)x_  eiK(n-2)xjGT>B, 

/ A  ^  1  C\ 


(4.2.15) 


where 


qT  =  pR^ _ dffi.  p 

Ba/T  Arr  O’ 


~  da 

(a)  (b)  (c) 


(4.2.16) 


Fo—  -i1f„ 


(b)  (c)  (d) 


(4.2.17) 


Figure  4.2  and  4.3  show  the  real  and  imaginary  parts  of  (4.2.16)  and  (4.2.17)  as  a 
function  of  wavenumber  for  U0  =  -0.5  (solid),  U0  =  0. 1  (long  dashed),  and  UD  =  0.5 
(short  dashed).  Figure  4.2  (a)  represents  G^,  (b)  represents  term  (b),  and  (c)  represents 
term  (c)  in  (4.2.16).  Figure  4.3  (a)  represents  GB,  (b)  represents  term  (b),  (c)  is  term  (c), 
and  (d)  is  term  (d)  in  (4.2. 17).  Comparison  of  Figures  4.2  and  4.3  shows  that  overall,  GB 
is  larger  in  magnitude  than  G^  and  the  contribution  from  each  of  the  terms  (b),  (c),  and  (d) 
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Wavenumber  (n)  Wavenumber  (n) 


Wavenumber  (n)  Wavenumber  (n) 


Figure  4.2  Plot  of  (4.2.16)  where  (a)  is  the  real  and  (b)  is  the  imaginary  contribution  to 
the  term  GT,  (c)  is  the  real  and  (d)  is  the  imaginary  contribution  to  GT  from  FB(9F0/9a), 

and  (e)  is  the  real  and  (0  is  the  imaginary  contribution  to  GT  from  F0(9Fg/9a)  as  a 
function  of  wavenumber.  The  solid  line  is  for  U0  =  -0.5,  long  dashed  is  for  U0  =  0.1,  and 
the  short  dashed  is  for  U0  =  0.5. 
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Figure  4.3  Plot  of  (4.2. 17)  where  (a)  is  the  real  and  (b)  is  the  imaginary  contribution  to 
the  term  G®,  (c)  is  the  real  and  (d)  is  the  imaginary  contribution  to  G®  from  Fg(5F0/9a), 

(e)  is  the  real  and  (0  is  the  imaginary  contribution  to  G®  from  FQ(9F3/9a),  and  (g)  is  the 

real  and  (h)  is  the  imaginary  contribution  to  G®  from  F0  as  a  function  of  wavenumber. 

The  solid  line  is  for  U0  =  -0.5,  long  dashed  is  for  U0  =  0. 1,  and  the  short  dashed  is  for  U0 
=  0.1. 
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is  larger  at  the  bottom.  Figure  4.3  (g)  and  (h)  shows  that  the  contribution  due  to  the 
perturbation  flowing  over  the  mountain  is  small  compared  the  contributions  from  the  other 
terms. 

The  two  terms  on  the  right  hand  side  of  (4.2.15)  involve  the  three  wavenumbers  n, 
n+2,  and  n-2  which  represent  the  forcing  of  0(v  *)  by  0(v°).  Therefore,  we  expect  X\ to 
be  forced  by  three  terms  proportional  to  these  wavenumbers.  We  will  assume  solutions  for 
X  \  of  the  form 

XI  =  Xl,neiKnx  +  Xl,n+2e'K<n+2>x  +  Xl,n-2eiK(n-2)x  +  CC. 

(4.2.18) 

Substituting  (4.2.18)  into  (4.2.15)  and  projecting  the  equation  onto  the  modes  n,  n+2,  and 
n-2,  we  obtain  three  sets  of  equations  for  0(v^).  The  equations  are 
for  mode  n 


^--s((Knp-f  /2)1(X1,„)  =0, 

daz 

[-(%  +  uT.BKnlj^)  -KnXxi„  =  co,^. 


(4.2.19) 

(4.2.20) 


and  for  mode  n+2 


(4.2.21) 


and  for  mode  n-2 


~  s((K(n+2))2+  fl)  |  (XI, n+2)  -0, 

[-do  +  UT-BK(n+2)] (8Xl^'1'2)  -XK(n+2)xl,n+2  +UKGT3.0, 

(4.2.22) 


•  s((K(n-2))2+  P) 

9<r 


{XI, n-2)  =0, 


(4.2.23) 
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[-<Bo  +  uT.BK(n-2)]  -XK(n-2)xi,n-2  +UKGT3*0. 

(4.2.24) 

At  0(v  1 ),  co  j  enters  the  problem  in  (4.2.20)  only.  To  evaluate  co  j ,  we  combine 
(4.2.20),  (4.2.19)  and  the  0(v°)  thermodynamic  boundary  condition,  (4.2.1 1)  and  its 
solution.  We  find  to  j  must  equal  zero  which  implies  that  by  proper  gauging,  we  can 
choose  Xin  =  0.  Asa  consequence  of  to  j  =  0,  the  0(v  ^)  problem  will  result  in  a 
correction  to  the  modal  structure  for  the  O(v^)  problem.  To  get  a  correction  to  the 
frequency,  co0,  we  must  solve  the  O(v^)  problem  and  find  0)2-  Therefore,  we  cannot  infer 
anything  about  the  group  velocity  from  these  two  orders  of  the  problem. 

We  have  separated  the  problem  into  0(v°)  and  O(v^).  After  solving  the  leading  order 
Eady  problem,  we  could  choose  some  wavenumber  n.  That  wavenumber  might  be  in  the 
neutral  or  unstable  region  shown  in  Figure  4.1.  If  n  is  in  the  unstable  region,  it  will  have 
growth  rate  imaginary(co0)  shown  in  Figure  4.1  (a)  and  propagate  at  frequency  real(to0) 
shown  in  Figure  4. 1  (b).  If  n  is  in  the  stable  region,  there  will  be  two  waves  that  move  at 
frequency  (0o+’*  shown  in  Figure  4.1  (b)  and  (c)  depending  on  the  sign  we  choose  in 
front  of  the  radical  in  (4.1.10). 

We  found  that  only  o>0  enters  the  0(v  * )  correction  to  the  Eady  problem.  Therefore,  the 
modes  n+2  and  n-2  will  propagate  and  grow  as  n.  If  the  wavenumber  n  is  in  the  unstable 
region  shown  in  Figure  4.1,  n+2  and  n-2  will  propagate  at  real(co0)  and  have  growth  rate 
imaginary(<0o).  If  the  wavenumber  n  is  in  the  stable  region,  n+2  and  n-2  will  propagate 
with  frequency  co0+»*  depending  on  the  sign  we  choose  in  (4.1.10). 

\Yc  solved  the  0(v  * )  problem  following  the  same  procedure  as  in  Section  4.1  for  the 
Eady  model.  Equations  (4.2.19),  (4.2.21),  and  (4.2.23)  have  solutions  of  the  form 

Xn,n±2<°> =  An,n±2cosh(  2^n,n±2(° '  2*) +  Bn,n±2sinh(  2^n,n±2<° '  §•>)  • 

(4.2.25) 
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The  thermodynamic  equations  at  the  upper  and  lower  boundaries,  (4.2.20),  (4.2.22),  and 
(4.2.24)  were  applied.  Since  (Dj  =  0,  we  can  solve  for  the  coefficients  An  and  Bn 
immediately  because  (4.2.19)  and  (4.2.20)  represent  the  Eady  model  which  we  already 
solved  in  Section  4.1.  Next,  we  wrote  a  system  of  algebraic  equations  for  the  unknown 
coefficients  An+2  and  Bn+2  and  solved  the  system  using  a  numerical  technique.  After  we 
found  the  coefficients,  we  reconstructed  %  i  >rl+2  and  finally  the  full  streamfunction. 

Figure  4.4  shows  the  amplitudes  of  the  coefficients  X0>  X\  n+2>  and  X\  n-2’ at  the  (a) 
lower  and  (b)  upper  boundaries.  For  this  case,  we  chose  the  Burger  number,  S  =  0.25  and 
the  most  unstable  mode  is  wavenumber  six.  At  the  0(v  ^),  the  modes  four  and  eight  enter 
the  problem.  The  magnitude  of  the  coefficients  are  larger  at  the  lower  boundary.  We  will 
compare  this  figure  to  results  from  the  numerical  solution  in  Section  4.3. 

Figure  4.5  shows  x-o  contours  of  the  perturbation  streamfunction  down  the  center  of 
the  channel.  Figure  4.5  (a)  shows  the  full  streamfunction  which  is  the  same  of  the  leading 
order  streamfunction  plus  the  correction  due  to  topography.  Figure  4.5  (b)  is  the 
streamfunction  for  the  leading  order  Eady  problem.  We  chose  the  stratification,  S  =  0.25 
and  the  fastest  growing  mode  is  wavenumber  six.  To  find  the  correction  to  the  Eady 
model,  we  chose  v  =  0.25.  Since  S  =  0.25,  the  ratio  of  mountain  height  to  Rossby 
number  is  one.  These  are  good  choices  to  remain  within  the  quasi-geostrophic 
approximation.  Figure  4.5  (c)  and  (d)  show  the  correction  due  to  wavenumber  eight  and 
wave  number  four.  All  three  modes  have  the  same  growth  rate  and  frequency. 

4.3  Numerical  Solution 


In  this  section  we  use  a  numerical  technique  to  solve  the  topography  problem.  First, 
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Wavenumber  (n) 


Figure  4.4  The  magnitude  of  the  coefficients  in  (4.2.14)  and  (4.2.18)  for  (a)  the  lower 
and  (b)  the  upper  boundary  for  the  case  of  n  =  6,  n  +  2  =  8,  and  n  -  2  =  4.  Here,  U„  = 
0.1,  e  =  0.1,  and  \  =  1.0. 


Figure  4.5  X-c  cross  sections  at  mid-radius  for  the  perturbation  streamfunction  for  (a) 
the  sum  of  the  most  unstable  mode  from  the  Eady  model,  wavenumber  6,  plus  corrections 
due  to  topography  from  wavenumbers  4  and  8.  (b)  Most  unstable  Eady  mode, 
wavenumber  six,  (c)  wavenumber  8,  and  (d)  wavenumber  4.  S  =  0.25,  U0=  0.0,  and  X 
=1.0.  Topographic  ridges  at  0.0  and  n  rad. 
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we  cast  the  equations  in  spectral  form.  Next,  we  manipulate  the  equations  to  set  up  an 
eigenvalue,  eigenvector  problem  and  solve  the  problem  using  numerical  routines.  The 
eigenvalues  give  us  the  growth  rates  and  frequencies  and  the  eigenvector  will  be  used  to 
reconstruct  the  streamfunction.  Finally,  we  consider  the  topographic  instability  problem. 

In  Section  4.2,  we  showed  the  quasi-geostrophic  potential  vorticity  equation  in  sigma 
coordinates  with  a  wavy  basic  state  to  be 


|UjOF,q)  =  0, 

with  the  boundary  condition  at  a  =  0.5 


djd(p 

at 


(4.2.1) 


(4.2.2) 


and  the  boundary  condition  at  c  =  -0.5 


8148145 


+  vr|  =  0, 


where 


'F(x,y,o,t)  =  -y(u0  +  X(o+  hj  +  vFB(o)cos(2Kx), 


(4.2.3) 


(4.3.1) 


and  v  =  (S8/e).  In  the  model,  there  is  a  rigid  upper  lid  and  cos(2Kx)  bottom  topography. 
We  take  the  mode!  to  be  periodic  in  the  x  and  y  directions.  We  chose  y  to  be  periodic  to 
avoid  the  difficult  problem  of  the  intersection  of  the  side  wall  with  the  topography.  In  this 
problem,  we  are  interested  in  the  phenomena  in  the  center  of  the  channel,  far  from  the 
walls.  Wc  assume  solutions  to  the  full  equations  of  the  form 

M  N 

<p(x,y,o,0  =  -u'(i)y  +  £  £  [<Pmn(i.a)eiKnx  +  <(taS(t,a)e'HsinMy)  + 

IT1=0  71=0 

M  N 

X  X  [<Pnin(t,a)eiKnx  +  cp^(t,o)e-lKnx]cos(m/y). 

711=0  71=0 


(4.3.2) 
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In  (4.3.2),  *  represents  the  complex  conjugate.  u'(t)  represents  a  volumetrically  averaged, 
time  dependent,  zonal  velocity.  It  becomes  important  when  we  consider  the  interaction  of 
the  zonal  flow  with  topography  through  the  form-drag  equation.  We  represent  periodicity 
in  y  using  a  sin(m/y)  and  cos(m/y)  rather  than  exp(im/y)  to  prevent  wave  propagation  in  the 
north-south  direction.  In  (4.3.2)  the  amplitudes  are  functions  of  time  and  the  vertical 
coordinate. 

As  in  the  Eady  model,  if  we  take  q  =  0  at  time  t  =  0,  then  q  remains  zero  for  all  time  and 
(4.2.1)  becomes 


(4.3.3) 


Substituting  (4.3.2)  into  (4.3.3),  the  discrete  solutions  are  of  the  form 


9mn(<*)  —  A^n(t)cosh(  2jImno)  +  Bfnn(t)sinh(  2jimno)  +  CC 


and 


(4.3.4) 


where 


<Pmn(°)  -  Arnn(t)cosh(  ^mn°) +  ®mn(0sinh(  2(Irano)  +  CC 

(4.3.5) 


^mn  =  ^(s((Kn)2  +  (/m)2))]. 


(4.3.6) 


Next,  we  substitute  the  assumed  solution,  (4.3.2),  into  the  thermodynamic  boundary 
conditions  (4.2.2)  and  (4.2.3)  and  obtain  expressions  involving  summations  of  zonal  and 
meridional  wavenumbers.  To  separate  the  equations  into  individual  zonal  and  meridional 
components,  we  project  the  equations  onto  orthogonal  basis  functions.  We  take  the  sine 
and  cosine  as  the  meridional  basis  functions  and  project  on  an  arbitrary  wavenumber  r.  We 
take  an  exponential  function  as  the  zonal  basis  function  and  project  on  an  arbitrary 
wavenumber  q.  Projecting  on  the  sine  and  cosine  in  the  meridional  direction  will  produce 


two  sets  of  equations.  After  some  algebraic  manipulation,  the  boundary  conditions  in 
spectral  form  for  the  cosine  projection  may  be  written 


(4.3.7) 

at  o  =  0.5  and 


(4.3.8) 

at  o  =  -0.5.  A  similar  set  of  equations  can  be  written  for  the  sine  component.  Here,  the 
superscripts  +  and  -  indicate  the  quantities  should  be  evaluated  at  the  upper  or  lower 
boundary.  There  are  two  complex  conjugate  pairs  with  the  superscript  s  in  (4.3.7)  and 

(4.3.8).  The  complex  conjugate  becomes  important  when  we  consider  zonal  wavenumbers 
zero,  one,  and  two. 

There  are  many  interesting  points  to  be  made  about  these  equations  and  the  same  things 
can  be  said  for  the  sine  component.  The  equations  are  decoupled  in  the  meridional 
direction  which  means  an  independent  set  of  equations  can  be  written  for  each  meridional 
wavenumber.  The  equations  are  coupled  in  the  zonal  direction  by  q±2  due  to  interaction  of 
the  perturbation  with  the  wavenumber  two  of  the  topography  and  the  forced  basic  state. 

The  equations  for  the  cosine  projection  involve  sine  and  cosine  coefficients.  The  sine 
coefficients  enter  because  we  took  y-derivatives  in  the  original  equation.  The  second  and 
fifth  terms  in  (4.3.7)  and  (4.3.8)  have  the  sine  coefficients.  The  second  term  in  both 
equations  represents  the  basic  state  advection  of  perturbation  temperature  with  wavenumber 
q±2.  The  fifth  term  in  both  equations  represents  advection  by  perturbation  wavenumber 
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q±2  of  the  basic  state  temperature.  Also,  the  fifth  term  in  the  lower  boundary  condition, 

(4.3.8),  represents  perturbation  advection  of  wavenumber  q±2  along  the  slope  of  the 
topography.  The  magnitudes  of  each  of  these  terms  is  controlled  by  the  ratio  v  =  (S8/e) . 
Thus  the  stratification,  mountain  height,  and  rotational  effects  will  all  play  a  role  in 
determining  how  important  the  wavy  basic  state  and  the  mountain  are  in  this  problem.  The 
third  and  fourth  terms  in  (4.3.7)  and  (4.3.8)  represent  the  Eady  model.  In  the  absence  of 
topography  and  wavy  basic  state,  this  model  reduces  to  the  Eady  model. 

To  solve  these  equations,  we  substitute  for  tpqr  in  (4.3.7)  and  (4.3.8)  using  (4.3.4)  and 
(4.3.5).  Next  we  add  and  subtract  the  upper  and  lower  boundary  conditions.  This  step 
allows  us  to  solve  for  the  time  derivatives  in  (4.3.7)  and  (4.3.8).  Finally,  after 
considerable  algebra,  we  derive  the  following  dynamical  system  of  equations  for  the 
unknown  coefficients 

Kf  +  i  [a:l(A2-qj  +  A2^J>r)  +  ac2A^+2,r  +  <as5Aq,r  +  6r5Bqj  +  be  l(B2^j  +  62^}+  bc2Bq+2j]  =  0. 

Bty  +  i  [cc3( Aq_2j  +  +  0C4A^j+2,r  +  ar6Aqj  +  bs6Bqj  +  bc^Bq_2j  +  B2-qj)+  bc4Bq+2j\  =  0, 

Aq-  +  i  [as1^  A^_2j  +  A2Hqj)  +  as2A^+2,r  +  uc’5Aqj-  +  bcSBqj  +  bs  l(Bq_2j  +  B2^j}+  £>r2Bq+2  j]  =  0, 

Btf  +  i  [av3(Acj_2j  +  A2^j)  +  as4A^2,r  +ac6Aqj  +  £>c6Bqj  +  bs^B^j  +  B2Sqj}+  Z?r4Bq+2j]  =  0. 

(4.3.9) 


All  of  the  coefficients  written  in  italics  are  listed  in  Appendix  B. 

The  system  represents  a  set  of  complex,  linear,  ordinary  differential  equations  and  can 
be  written  as 


(4.3.10) 


where  ©  is  a  vector  containing  the  Fourier  coefficients  and  A  is  a  square  coefficient  matrix. 
Since  this  is  a  linear  system,  we  can  assume  exp(cot)  time  dependence  and  transform  the 
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system  to  the  eigenvalue,  eigenvector  problem, 

to©  =  A©.  (4.3.11) 

The  eigenvalue,  to,  and  the  eigenvector  are  computed  using  IMSL  library  routines.  In  this 
problem,  real(to)  represents  the  growth  rate  and  imaginary(to)  represents  the  frequency. 
Associated  with  each  eigenvalue  is  an  eigenvector  that  will  be  used  to  reconstruct  the 
streamfunction  using  an  inverse  Fourier  transform  routine. 

We  truncated  the  system  at  zonal  wavenumber  twelve  and  took  the  meridional 
wavenumber  to  be  one.  We  expect  the  largest  growth  rates  at  approximately  wavenumber 
six  and  since  wave  six  interacts  with  waves  eight  and  four,  we  should  include  at  least 
through  wavenumber  eight.  However,  Niehaus  (1980)  showed  that  growth  rates 
continued  to  change  as  the  truncation  level  was  increased  until  twelve  or  more  waves  were 
added  to  the  model.  Wavenumbers  greater  than  twelve  were  neglected.  Wavenumbers 
zero,  one,  and  two  were  treated  as  special  cases.  Wavenumber  zero,  the  zonal  mean,  and 
one  interact  with  the  complex  conjugate  of  wavenumbers  two  and  one.  Wavenumber  two 
interacts  with  the  mean  zonal  flow. 

Figure  4.6  (a)  is  a  plot  of  contours  of  growth  rate  for  the  fastest  growing  mode  and  (b) 
is  a  plot  of  contours  of  frequency  as  the  stratification  and  mountain  height  are  varied. 

When  the  mountain  height  is  zero,  the  ordinate  gives  the  growth  rates  of  the  Eady  model, 
and  Figure  4.6  (a)  shows  that  as  the  stratification  increases,  the  growth  rate  decreases.  As 
the  mountain  height  increases,  the  growth  rates  decrease  slightly  for  large  values  of 
stratification  but  are  nearly  constant  for  small  values  of  stratification.  Figure  4.6  (b)  shows 
that  the  frequency  remains  nearly  constant  until  larger  mountain  heights  and  then  the 
frequency  begins  to  increase.  In  Figure  4.6,  £  =  0.1,  UQ  =  0. 1,  and  X  =  1.0.  When  UQ 
=  -0. 1,  a  plot  similar  to  Figure  4.6  shows  that  the  growth  rate  and  frequency  remain  nearly 


Figure  4.6  Contours  of  (a)  growth  rate  and  (b)  frequency  for  the  most  unstable  mode 
a  function  stratification,  S,  and  mountain  height,  li.  Allows  indicate  the  direction  of 
increasing  values  of  growth  rate  or  frequency.  U0=  0. 1 ,  X  =  1 .0,  and  £  =  0.1. 
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constant  for  all  mountain  heights  and  stratification.  For  all  of  the  remaining  figures,  the 
vertical  shear,  X  -  1.0. 

Figure  4.7  (a)  is  a  plot  of  contours  of  growth  rate  for  the  fastest  growing  mode  and  (b) 
is  contours  of  frequency  as  the  stratification  and  U0  are  varied.  UQ  varies  on  the  abscissa. 
The  mountain  height  and  Rossby  number  are  fixed  at  0. 1.  Figure  4.7  (a)  shows  the 
growth  rate  of  the  most  unstable  modes  is  nearly  constant  for  U0  >  -0.1.  As  UQ 
approaches  -0.5,  the  growth  rate  becomes  extremely  large.  Figure  4.7  (b)  shows  the 
frequency  increases  as  UQ  increases  but  decreases  as  the  stratification  increases.  As  UQ 
approaches  -0.5,  the  frequency  of  the  mode  becomes  zero.  This  happens  since  the  phase 
speed  of  propagation  is  equal  and  opposite  to  U0.  Although  this  mode  is  growing  but  has 
no  frequency,  it  should  not  be  considered  a  form  of  topographic  instability  since  the  form 
drag  mechanism  was  not  considered.  Surface  easterly  flows  of  such  magnitude  are  not 
observed  in  the  atmosphere  or  the  annulus  so  this  flow  is  probably  not  physically 
realizable.  As  the  mountain  height  was  increased,  plots  similar  to  Figure  4.6  showed 
nearly  the  same  growth  rates  for  positive  U0  but  growth  rate  becomes  enormous  for  all 
negative  values  of  UQ. 

We  wanted  to  investigate  the  effect  of  the  mountain  height  on  the  basic  state  and 
perturbation.  We  fixed  S  =  0.25  and  e  =  0.1.  Figure  4.8  shows  x-y  and  x-a  cross  sections 
of  the  basic  state  and  perturbation  stream  function.  Figure  4.8  (a)  and  (b)  show  that  for 
very  small  mountain  heights,  the  basic  state  flow  at  the  bottom  and  top  is  nearly  zonal. 
When  the  x-dependent  part  is  added  to  the  zonal  flow,  only  a  linear  vertical  shear  is 
apparent  in  Figure  4.8  (c).  Figure  4.8  (d)  is  the  x-a  cross  section  of  the  x-dependent  part 
of  the  basic  state  shows  a  high  over  the  topography  at  0  and  K  rad.  Figure  4.8  (e)  and  (0 
show  the  perturbation  streamfunction  at  the  bottom  and  top  for  the  same  small  mountain. 
The  perturbations  line  up  zonally  much  as  is  observed  in  an  Eady  model  which  is  periodic 
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Figure  4.7  Contours  of  (a)  growth  rate  and  (b)  frequency  for  the  most  unstable  mode  as 
a  function  of  S  and  UQ.  Arrow  indicates  the  direction  of  increasing  values.  8  =  0.1,  A.  = 
1.0,  and  e  =  0.1. 
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Figure  4.8  Basic  state  and  perturbation  streamfunction  for  very  small  mountain  height,  h 
=  0.01.  (a)  x-y  cross  section  of  basic  state  streamfunction  at  top  and  (b)  bottom,  (c)  x-c 
cross  section  of  basic  state  streamfunction  at  mid-radius  and  (d)  wavy  part  of  basic  state 
streamfunction.  (e)  x-y  cross  section  of  perturbation  streamfunction  at  top  and  (b)  bottom, 
(g)  x-a  cross  section  perturbation  at  mid-radius.  Topographic  ridge  at  0.0  and  rc  rad.  Uq 
=  0. 1 ,  X  =  1.0,  and  e  =  0. 1 . 


75 


in  y.  Figure  4.8  (g)  is  an  x-c  cross  section  of  the  perturbation  and  is  nearly  the  same  as  the 
perturbation  for  the  Eady  model  shown  in  Figure  4.5  (b). 

Next,  we  increased  the  mountain  height  to  0.1.  Figure  3.3  (a)  and  (b)  shows  the  basic 
state  streamfunction  is  very  wavy  at  the  bottom  and  top.  Figure  3.3  (c)  show  the  x-ct  cross 
section  is  also  wavy  with  a  high  pressure  ridge  over  the  mountain  ridge  at  0  and  k  rad. 
Figure  4.9  (a)  and  (b)  shows  the  x-y  cross  section  of  the  perturbation  streamfunction  at  the 
top  and  bottom.  The  perturbations  follow  the  wavenumber  two  basic  state.  Thus,  the 
basic  state  streamfunction  acts  as  a  guide  or  ‘storm  track’  for  the  perturbation.  Figure  4.9 
(c)  shows  an  x-a  cross  section  of  the  perturbation  streamfunction.  As  predicted  from  the 
analytical  model,  there  is  one  dominant  wave  number  but  the  width  of  the  highs  and  lows 
shows  there  is  a  superposition  of  several  waves.  Just  downstream  from  the  mountain 
ridge,  the  upper  and  lower  modes  are  closer  together  and  ‘feel’  each  others  presence  and 
they  have  large  amplitude  at  the  top  and  bottom.  Thus,  we  expect  enhanced  cyclogenesis  in 
this  region.  Figure  2.5  (b)  shows  that  downstream  of  the  Rocky  Mountains,  enhanced 
cyclogenesis  occurs. 

Figure  4.10  shows  the  amplitudes  of  the  spectral  coefficients  of  (4.3.2)  for  the  most 
unstable  mode  for  the  same  parameters  values  considered  in  figure  4.9  as  well  as  in  Figure 
4.4  for  the  analytical  solution.  Wavenumbers  four,  six,  and  eight  have  magnitudes  about 
an  order  of  magnitude  larger  than  the  other  wavenumbers  for  this  mode.  These  are  also  the 
most  important  wavenumbers  from  the  analytical  model.  As  in  Figure  4.4, 4.10  shows  the 
magnitudes  of  the  coefficients  at  the  lower  boundary  have  the  largest  amplitudes. 

We  also  investigated  the  effect  that  U0  had  on  the  basic  state  and  perturbations.  When 
U0  is  negative,  the  basic  state  is  out  of  phase  with  the  mountain.  This  is  not  observed  in 
the  atmosphere  or  the  annulus  and  probably  is  not  physically  realizable.  When  UD  is 
positive,  the  basic  state  is  in  phase  with  the  mountain  and  shown  Figure  3.3  (a)  -  (d). 
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Figure  4.10  Magnitude  of  the  coefficients  of  (4.3.2)  as  a  function  of  wavenumber  for 
(a)  the  lower  and  (b)  the  upper  boundary.  Here,  S  =  0.25,  Un=  0.1,  8=  0.1,  e  =  0.1,  and 

A.  =  1.0. 
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When  UQ  =  0.0,  there  is  no  x-dependence  in  the  basic  state  but  the  perturbation  can  still 
interact  with  the  mountain.  Figure  4.1 1  (a)  and  (b)  shows  x-y  cross  sections  at  the  top  and 
bottom.  As  the  perturbation  flows  over  the  mountain  it  moves  to  the  north.  This  can  be 
understood  form  the  thermodynamic  boundary  condition,  (4.2.1 1).  As  the  zonal  flow  rises 
over  the  mountain  slope,  it  cools  adiabatically.  The  response  in  the  fluid  is  an  increased 
meridional  flow  which  advects  warmer  fluid  and  the  perturbations  north.  Figure  4.11  (c) 
shows  the  vertical  structure  of  the  perturbation.  Buzzi  (1987)  showed  the  same  type  of 
phenomena  occurs  in  their  two-layer  quasi- geostrophic  model.  Buzzi  (1989)  showed  that 
as  observed  perturbations  flow  over  the  mountain,  they  move  north  and  rotate 
anticyclonically  as  shown  in  Figure  2.6. 

To  investigate  the  time  dependence  in  the  continuous  model,  we  plotted  the  x-z  and  x-y 
perturbation  contours  as  a  function  of  time.  The  plots  shown  in  Figure  4.12  (a)  -  (d)  are  x- 
a  cross  sections  and  (e)  -  (f)  are  x  -y  cross  sections  for  non-dimensional  times  of 
approximately  0.5  days.  UQ  =  0.0  so  there  is  no  wavy  basic  state.  At  each  time  step  we 
normalized  by  the  exponential  growth  rate.  In  Figure  4.12  (a)  -  (d),  the  x-a  cross  sections 
show  that  as  the  perturbation  rises  up  the  mountain  slope,  it’s  amplitude  is  decreased  at  the 
surface.  As  the  perturbation  flows  down  the  mountain  slope,  the  amplitude  at  the  surface  is 
enhanced.  Just  downstream  from  the  mountain,  instabilities  at  the  top  and  bottom  extend 
toward  each  other  suggesting  the  possibility  for  enhanced  growth.  Figure  2.5  taken  form 
observations  shows  this  is  a  region  of  enhanced  growth.  Figure  4.12  (e)  -  (g)  shows  the 
x-y  cross  section  as  a  function  of  time.  As  anticipated,  the  perturbation  moves  north  as  it 
flows  over  the  mountain. 

In  Section  4.2,  the  analytical  model  showed  that  to  leading  order,  the  problem  is  the 
Eady  model  with  a  correction  due  to  the  mountain.  Therefore,  we  should  expect  in  this 
numerical  model,  that  the  fastest  growing  mode  will  be  that  predicted  by  the  Eady  model 


Figure  4.11  Perturbation  streamfunction  when  the  basic  flow  is  zero  so  that  there  is  no 
wavy  basic  state  but  the  mountain  height  is  h  =  0.1.  (a)  x-y  cross  section  at  top  and  (b) 
bottom,  (c)  x-o  cross  section  at  mid-radius.  Topographic  ridge  at  0.0  and  n  rad.  S  = 
0.25,  Un  =  0.0.  e  =  0. 1 ,  and  1.0. 


Figure  4.12  Perturbation  streamfunction  as  a  function  of  non-dimensional  time  with  no 
wavy  basic  state,  (a)  -  (d)  x  a  cross  sections  and  (e)  -  (0  x-y  cross  sections.  Topographic 
ridge  at  0.0  and  n  rad.  S  =  0.25,  U0  =  0.0,  e  =  0.1,  and  X=  1.0. 


Figure  4.12  (continued) 
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and  if  S  =  0.25,  the  wavenumber  is  six.  The  mode  with  the  second  largest  growth  could 
be  found  using  Figure  4.1.  and  this  mode  matches  the  mode  found  in  the  numerical  model, 
The  frequency  for  every  mode  found  in  the  numerical  model  matches  that  found  with  the 
Eady  model.  Therefore,  the  analytical  model  provides  a  way  to  understand  the  numerical 
results. 


When  we  write  the  full  set  of  equations  in  spectral  form,  we  find  that  the  zonal 


acceleration  of  the  mean  flow  enters  the  equations  when  the  meridional  wavenumber  is 
zero.  When  the  meridional  structure  is  zero,  we  derive  a  different  set  of  spectral  equations 


(U0+XXi2K)^-  -  (Xi2K)<p++  ivKvF+  u\2  =0, 


(4.4.1) 


|fe]  +  uWg-  -  (X2kW+  ivK('P'+  l)  uA,,2  =  0. 


(4.4.2) 


In  (4.4.1)  and  (4.4.2),  the  last  term  represents  the  advection  of  the  volumetrically 


averaged  u'(t)  by  the  basic  state  and  over  the  mountain  slope  at  the  lower  boundary  and  8 


represents  the  Kronecker  delta.  This  term  appears  only  when  the  perturbation  has  the  same 
wavenumber  as  the  basic  state  or  mountain.  Therefore,  these  equations  apply  for  zonal 


wavenumber  two  only.  For  other  perturbation  wavenumbers,  the  last  term  is  zero  and  the 
problem  reduces  to  the  Eady  problem.  By  introducing  a  new  variable  u'(t)  into  these 
equations,  we  will  need  one  additional  equation,  the  form-drag  equation,  to  solve  the 
problem.  If  we  substitute  the  wavy  basic  state  (4.3.1)  and  perturbation  (4.3.2)  into  the 
form-drag  equation,  (3.7.4),  after  some  algebra  we  obtain  an  expression 


83 


«<t)  =  ^(|)f  AUcosh(ji20) '  Sin~~2Q^|  *  BI20sinh(|i20) . 

®Ve/L  \  ^20  /  J  (4.4.3) 

^20  and  ®20  are  the  imaginary  parts  of  complex  coefficients.  For  the  problem  with  no  y- 
s tincture  we  have  three  complex  equations  and  three  unknowns. 

Buzzi  etal.  (1984,  hereafter  BTS)  also  solved  this  problem  with  no  y-structure  using  a 
continuous  Eady  model  and  wavy  and  zonal  basic  states.  However  in  their  model  as  they 
varied  the  wavenumber,  the  wavenumber  of  the  basic  state,  mountain,  and  perturbation  all 
changed.  We  cannot  make  a  direct  comparison  of  our  work  with  the  work  of  BTS  since 
our  basic  state  and  mountain  are  always  wavenumber  two  but  we  can  compare  when  the 
perturbations  is  wavenumber  two.  Also,  our  model  does  not  include  friction. 

To  solve  (4.4.1)  and  (4.4.2)  we  follow  the  same  procedure  outlined  in  Section  4.2. 
Here  we  only  consider  wavenumber  two,  the  wavenumber  of  the  topography  and  the  basic 
state.  When  the  form-drag  is  included,  the  system  of  equations  is 


B  +  alA  +  MB  +  ulu  =  0 

(4.4.4) 

A  +  a2A  +  b2B  +  u2u  =  0 

(4.4.5) 

u  +  a3A  +  b3B  =  0 

(4.4.6) 

The  coefficients  marked  with  italics  are  listed  in  Appendix  B. 

Figure  4.13  (a)  shows  contours  of  growth  rate  and  (b)  frequency  as  a  function  of 
stratification  and  UD.  We  find  topographic  instability  for  small  values  of  stratification  and 
for  positive  UQ.  Figure  4. 14  (a)  shows  contours  of  growth  rate  and  (b)  frequency  as  a 
function  of  stratification  and  mountain  height.  Figure  4.14  (b)  shows  that  as  the  mountain 
height  increases,  the  region  of  topographic  instability  occurs  for  larger  values  of 
stratification. 
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Figure  4.14  Contours  of  (a)  growth  rate  and  (b)  frequency  for  the  topographic 
instability  problem  as  a  function  of  S  and  h.  Region  of  zero  frequency  is  for  low  values  of 
S.  U0=  0.1,  X=  1.0  and  E  =  0.1. 


CHAPTER  5 


DISCUSSION  AND  CONCLUSIONS 

In  this  thesis,  we  developed  a  mathematical  model  to  describe  baroclinic  flow  over 
topography  on  a  wavy  basic  state.  We  chose  a  wavy  basic  state  based  on  the  observed, 
nearly  stationary  pattern  in  the  atmosphere.  We  were  more  interested  in  understanding  how 
baroclinic  waves  were  modified  by  topography  rather  than  topographic  instability  which  at 
present  appears  not  to  be  important  in  the  primitive  equation  model  of  Revell  and  Hoskins 
(1984). 

We  transformed  the  five  governing  equations  to  sigma  coordinates  so  that  the  lower 
boundary  coincided  with  the  lowest  sigma  level  in  the  model.  Based  on  scaling  arguments, 
we  derived  the  quasi-geostrophic  potential  vorticity  equation  and  thermodynamic  equations 
for  the  upper  and  lower  boundaries  in  sigma  coordinates.  In  sigma  coordinates,  the  effect 
of  topography  appears  in  the  thermodynamic  equation  at  the  lower  boundary.  The 
topography  term  is  proportional  to  v  =  S6/e,  or  the  mountain  height,  stratification  and 
Rossby.  We  purposely  chose  a  simple  model  with  only  the  essential  features  required  to 
describe  the  interaction  of  the  flow  with  the  mountain. 

To  solve  the  quasi-geostrophic  potential  vorticity  equation  and  the  thermodynamic 
equations  on  the  boundaries,  we  linearized  about  a  steady,  wavy  basic  state  with  linear 
vertical  shear  which  is  an  exact  solution  to  the  non-linear  problem.  The  x-dependent  part  of 
the  basic  state  was  forced  by  topography  and  simulates  the  observed  wavenumber  two 
standing  pattern  in  the  atmosphere  and  the  annulus.  We  have  made  the  simplification  that 
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both  the  basic  state  and  the  perturbation  have  zero  potential  vorticity  in  the  interior  of  the 
fluid,  in  the  spirit  of  the  Eady  model.  Nevertheless,  baroclinic  instability  is  possible  due  to 
the  potential  vorticity  sheets  found  on  the  horizontal  boundaries. 

We  solved  the  linear  instability  problem  using  an  analytical  and  a  numerical  technique. 
The  analytical  solution  capitalized  on  the  fact  that  all  the  terms  in  the  equation  involving  the 
mountain  or  the  wavy  basic  state  are  proportional  to  the  parameter  v  =  (S5/e)  which  was 
taken  to  be  small.  We  wrote  the  problem  using  an  asymptotic  expansion  in  that  small 
parameter  v  and  found  that  to  leading  order,  we  recover  the  Eady  model.  This  was  no 
surprise  since  to  0(v°),  the  basic  state  is  rectilinear  and  identical  to  the  Eady  basic  state  At 
the  0(v  1),  we  brought  in  the  correction  due  to  topography  and  wavy  basic  state.  A  very 
important  result  from  the  asymptotic  analysis  is  the  O(v^)  problem  provides  us  with  a 
correction  to  the  eigenfunction  structure  while  the  frequency  correction  is  only  found  at 
0(v^),  an  order  not  yet  considered.  We  selected  the  fastest  growing  mode  with 
wavenumber  n  to  correspond  to  the  leading  mode  n  of  the  Eady  model.  The  correction  to 
the  leading  order  problem  involved  modes  n±2  due  to  the  x-dependence  introduced  in  the 
lower  and  upper  boundary  conditions  by  the  topography  and  basic  state  which  are 
wavenumber  two.  All  of  the  modes  n,  n±2  have  the  same  frequency,  ojq. 

In  the  numerical  technique,  the  quasi-geostrophic  potential  vorticity  equation  and 
thermodynamic  boundary  equations  were  transformed  to  spectral  space.  Since  we  chose 
the  potential  vorticity  on  the  interior  of  the  fluid  equal  to  zero,  the  solution  in  the  vertical 
involved  hyperbolic  functions.  We  chose  the  model  to  be  periodic  in  x  and  y  so  that 
trigonometric  functions  were  solutions  in  the  horizontal.  The  transformation  resulted  in 
spectral  equations  that  were  coupled  in  the  zonal  direction  with  wavenumber  two  due  to 
topography  and  wavy  basic  state.  The  equations  were  decoupled  in  the  meridional 
direction.  We  separated  the  equations  into  a  dynamical  system  of  equations  for  the  time 
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dependent  spectral  coefficients.  Since  the  model  is  linear,  we  assumed  an  exponential  time 
dependence  and  solved  the  system  using  an  eigenanalysis.  The  analysis  returns  growth 
rates  and  frequencies  for  the  various  modes  in  the  problem  as  well  as  the  eigenvector  to 
reconstruct  the  streamfunction  for  a  given  mode. 

We  compared  the  frequencies  of  this  model  with  those  of  the  analytical,  Eady  model. 

We  showed  the  growth  rate  and  frequency  for  the  most  unstable  mode  found  in  the 
numerical  model  as  a  function  of  stratification  and  mountain  height.  If  the  mountain  height 
is  zero,  the  frequencies  are  those  of  the  Eady  model.  For  small  mountain  heights  and 
stratification,  v  =  S8/e,  is  small  and  the  analytical  model  and  the  numerical  results  are  in 
good  agreement.  However,  as  the  mountain  height  and  stratification  increase,  the 
frequency  predicted  by  the  analytical  model  departs  from  the  frequency  found  in  the 
numerical  calculation.  One  reason  for  the  divergence  is  the  analytical  model  was  based  on 
v  «  0(1)  and  another  is  there  were  many  more  modes  included  in  the  numerical  solution. 

The  growth  rate  was  found  to  be  sensitive  to  UQ  which  is  a  parameterization  of  the 
momentum  Lux  convergence  into  the  midlatitudes.  For  UQ  >-0.1,  the  growth  rate  is 
nearly  constant  as  UD  increases  and  the  frequency  appears  to  increase  linearly.  As  UQ  -> 
-0.5,  the  frequency  becomes  zero  and  the  growth  rates  become  enormous.  We  can 
understand  the  frequency  going  to  zero  from  the  Eady  dispersion  relation,  (4.1.10).  Since 
X./2  =  0.5,  the  frequency  approaches  zero  as  U0  approaches  -0.5.  In  Section  3.6,  we  show 
that  the  denominator  of  the  basic  state  coefficients  approaches  zero  as  UQ  approaches  -0.5 
which  causes  the  growth  rates  to  become  very  large. 

We  also  compared  what  happened  in  the  instability  problem  when  there  was  and  was 
not  a  wavy  basic  state.  When  U0  =  0,  the  wavy  part  of  the  basic  state  is  zero  and  there  is 
only  linear  vertical  shear  left  in  the  basic  state.  We  found  that  the  perturbations  tend  to 
follow  the  basic  state.  Thus,  the  basic  state  acts  as  a  ‘storm  track’  for  the  disturbances. 
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The  amplitude  of  the  basic  state  depends  on  UQ.  If  UQ  is  very  small  or  equal  to  zero  then 
the  basic  state  is  nearly  zonal.  But  the  perturbations  still  move  north  over  the  mountain. 
However,  if  U0  is  strong,  the  basic  state  has  large  north-south  excursions  and  the 
perturbations  are  carried  farther  north  than  they  would  have  been  without  a  wavy  basic 
state.  When  UQ  is  large,  more  of  the  basic  flow  is  being  forced  to  rise  over  the  mountain. 
As  it  rises,  it  cools  adiabatically  and  this  results  in  an  imbalance  in  the  thermodynamic 
equation.  To  compensate  for  the  imbalance,  a  southerly,  meridional  flow  is  established 
which  advects  warm  air  and  the  perturbations  from  the  south.  Observations  from  the 
atmosphere  show  that  the  baroclinic  eddies  propagate  north  as  they  flow  over  the  mountain. 
Also,  observations  from  the  atmosphere  show  the  storm  tracks  are  strong  on  either  side  of 
the  mountain  but  are  at  a  minimum  over  the  mountain.  Part  of  the  reason  for  this  may  be 
the  amplitude  of  the  basic  current  is  highly  dependent  on  UQ. 

Topographic  instability  occurs  in  this  model  when  there  is  no  y-structure  and  when  the 
perturbation  is  wavenumber  two,  the  wavenumber  of  the  basic  state  and  the  mountain.  It 
occurs  due  to  the  form-drag  which  couples  the  acceleration  of  the  mean  zonal  flow  with  the 
topography.  Revell  and  Hoskins  (1984)  questioned  whether  topographic  instability  exists 
in  a  purely  baroclinic  model.  Kallen  (1982)  and  Buzzi  et  al.  (1984)  showed  that  a 
baroclinic  model  with  constant  zonal  flow  does  have  the  mechanism  for  topographic 
instability  to  occur.  In  this  baroclinic  model,  the  constant  flow  UQ  is  present  and  Figures 
4.9  and  4. 10  show  the  region  in  parameter  space  where  it  occurs.  The  instability  occurs 
for  physically  unrealistic  values  of  stratification.  For  values  of  stratification  this  small,  the 
whole  atmosphere  would  be  unstable. 

We  chose  to  work  with  a  continuous  model  in  the  vertical  but  could  have  worked  with  a 
two  layer  model  as  many  modelers  have  done.  The  continuous  model  is  better  since  it 
provides  for  better  vertical  resolution.  Also,  the  conditions  for  instability  may  be  different 
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in  a  continuous  baroclinic  than  a  two-layer  baroclinic  model. 

The  most  serious  weakness  of  this  model  is  the  use  of  the  f-plane  rather  than  the  beta- 
plane.  The  f-plane  is  physically  unrealistic  for  the  planetary  scale  waves  we  considered  for 
the  basic  state.  A  significant  improvement  could  be  made  if  we  solved  the  Chamey  (1947) 
problem  rather  than  the  Eady  problem.  The  Chamey  model  would  allow  us  still  to  use  a 
continuous  model  in  the  vertical  but  on  a  beta-plane;  analytical  solutions  of  this  problem 
probably  would  involve  prohibitive  amounts  of  algebra  since  the  wavy  basic  state  will  now 
possess  potential  vorticity. 

If  the  beta  effect  and  friction  were  included  in  the  basic  state  calculation,  the  phase  of 
the  basic  state  coi>'d  be  adjusted  with  respect  to  the  mountain  ridge  rather  than  the  high 
pressure  being  coincident  with  the  mountain  ridge.  If  the  beta  effect  and  friction  were 
included  in  the  perturbation  problem,  there  could  be  a  long  wave  cutoff  due  to  the  beta 
effect  and  a  region  of  stability  added  due  to  friction. 

A  still  more  realistic  basic  state  would  have  horizontal  as  well  as  vertical  shear  but  this 
would  make  an  analytical  solution  more  difficult.  The  horizontal  shear  would  act  as  a  guide 
for  the  instabilities  much  as  the  basic  state  does  in  this  model.  However,  the  horizontal 
shear  eliminates  the  need  for  periodicity  in  y.  Side  wall  boundary  conditions  would  confine 
the  north-south  amplitude  of  the  wave  and  channel  the  flow  instabilities.  This  would  give 
something  more  like  observed  in  the  annulus  but  may  not  represent  the  atmosphere  since 
there  are  no  sidewalls  in  the  atmosphere  to  act  as  wave  guides. 

The  linear  model  we  used  allows  for  zonal  perturbation  wave  interactions  of  plus  or 
minus  two  due  to  the  topography  and  the  basic  state.  Also,  the  linear  model  shows  there  is 
an  acceleration  of  the  zonal  mean  flow  due  to  the  form-drag  mechanism.  However,  if  we 
included  non-linear  interactions  we  could  consider  wave-zonal  mean  flow  interactions  for 
all  wavenumbers.  Non-linear  interactions  could  help  to  account  for  the  actual  observed 


patterns  in  the  atmosphere  that  cannot  be  explained  by  linear  theory. 


APPENDIX  A 


SCALING  OF  CURVATURE  TERMS  AND  THE  BETA  EFFECT 

A.l  Curvature  Terms 


In  this  appendix,  we  want  to  show  the  relative  magnitude  of  the  curvature  terms  and  the 
beta-effect  to  other  terms  in  the  momentum  equations.  Pedlosky  (1987)  gives  die  non- 
dimensional  form  of  the  x-component  of  the  momentum  equation  in  spherical  coordinates 
as 


e  +  i={5uw  -  uv  tan(0)]l  -  v  --  +  8w-Cos^~-  = 
'dt  r*  '  sin(0o)  cos(0o) 


cos(0Q)  rQ  dp  I  1  | 

cos(0)  r*  dx  \l+EpFj  (a  ]) 


where 


d  =  d  +ucos(80)r0  a  r0  g_  a 

dt  A  cos(0)  r* dx  r*  dy  dz  ' 


F=(y£ 

gD  * 

0  =  latitude, 


0O  =  mid-latitude  reference, 


and 


rG  =  earth's  radius. 
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There  are  no  approximations  to  the  momentum  equation.  The  equations  are  scaled  so  that 
the  relative  magnitude  of  each  term  may  be  measured  by  the  term  multiplying  it.  The  first 
two  terms  on  the  left  hand  side  of  (A.l)  are  the  advection  and  curvature  terms.  Each  of 
these  terms  is  scaled  by  the  Rossby  number.  If  we  decide  to  keep  terms  of  the  0(e),  we 
must  decide  on  the  relative  magnitude  of  the  advection  compared  to  curvature.  This  is 
measured  by  L/r*. 

The  circumference  of  the  Earth  at  a  given  latitude  is  (2ttr*  tan(9)).  If  we  choose  mid¬ 
latitude  to  ?  45°N,  the  length  of  fhe  Earth  is  2rcr*.  We  are  considering  wavenumber  two. 
Therefore,  the  length  of  one  wave  is  rcr*.  The  length  scale  for  this  wave  is  r*/2. 
Therefore,  a  measure  of  die  magnitude  of  curvature  compared  to  advection  is 


Tnus,  the  curvature  term  is  half  the  magnitude  of  the  advective  term.  As  a  first 
approximation,  we  can  neglect  this  term. 

For  the  annulus.  Bird,  Stewart,  and  Lightfoot  (1960)  give  the  x-component  of  the 
horizontal  momentum  equation  in  cylindrical  coordinates  and  non-dimensional  form  as 


where  R  is  the  radius  of  the  annulus.  Again,  if  we  are  to  include  terms  on  the  0(e),  then 
we  need  to  compare  the  magnitude  of  the  advective  terms  to  the  curvature  terms.  In  Section 
4.3  we  find  the  circumference  of  the  annulus  at  mid-radius  is  7t(a  +  b).  The  length  of  one 
mountains  is  half  the  length  of  the  channel.  Thus,  L  =  (a  +  b)/4  =  5.55  cm  and  the  radius  of 
the  channel  is  (b  +  a)/2.  Thus  the  measure  of  the  magnitude  of  curvature  compared  to 
advective  terms  is 


L_  (a  +  b)  __  5 
R  4(a  +  b)/2  ' 
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Therefore,  to  a  first  approximation  we  can  neglect  the  curvature  term  over  the  advection 
term. 


A. 2  Beta  Effect 

In  the  atmosphere,  the  depth  of  the  fluid  is  small  compared  to  the  horizontal  length  scale 
so  that  the  normal  component  of  the  Earth’s  vorticity  is  dynamically  significant.  Also,  the 
north-south  scale  of  the  motion  will  be  considered  small  so  the  Earth  can  be  considered  flat. 
The  variation  of  the  coriolis  parameter  with  latitude  can  be  linearized  about  a  mean  latitude 
so  that  f  may  be  written  as 

f  =  fo  +  Poy-  (A. 3) 

This  is  the  beta  -  plane  approximation.  In  our  model  we  chose  to  use  an  f  -  plane  which 
fixes  the  rotational  effect  of  the  Earth  as  constant.  A  check  of  the  relative  magnitudes  of  the 
terms  in  (A.  3)  gives 

2flcos(80)  L 

= - * - =  ^an(0o) . 

fo  2Qsin(0o)  0  (A.4) 

If  0O=  45°N,  Ur0=  .5  for  the  atmosphere.  In  our  model,  we  included  terms  on  the  0(e). 

For  planetary  waves,  an  estimate  of  the  Rossby  number  is  given  by 

211  = _ 2*.  .10  =003 

(6.3  x  106)  *  10"4 

In  the  annulus  experiments  we  want  to  simulate,  the  bottom  is  not  radially  dependent  so 
there  is  no  beta-effect.  However,  there  is  a  simulated  beta-effect  since  the  upper  surface  of 
the  fluid  slopes.  The  beta-effect  can  be  estimated  by  making  a  balance  of  forces  on  a  fluid 
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parcel  on  the  sloping  upper  surface  where 

R_  f  dh  2fl3R 

P  D3R  gD  ' 

The  magnitude  of  this  term  is  less  than  0(e)  and  is  neglected. 


(A.5) 


APPENDIX  B 


COEFFICIENTS  FOR  THE  NUMERICAL  SOLUTION  AND 
TOPOGRAPHIC  INSTABILITY  PROBLEM 

EJL _ Numerical  Solution 


Coefficients  for  (4.3.9)  are  as  follows: 

=  { M^q-2 jr  Ab  cosh(|l2o)  -  [M-20  Ab  sinh(^o)  -  -5]  coth(Hq-2.r)}  ~  , 

sinntpqj;  Hqj 

as2  =  {-  Mq+2/  Ab  cosh(|A2o)  +  [M-20  Ab  sinhto)  -  -5]  coth(Mq+2j)}  ^  , 

bsl  =  {[iq-2j  BB  sinh(M20)  coth(Mq.2/)  -  (>20  Bb  cosh(n2o)  +  .5]  }  ^  > 

sinnt|jq^  (iqj 

bsl  =  {-  liq+2>r  Bb  sinh(n20)  coth((iq+2j)  +  [M-20  Bb  cosh(|i2o)  +  .5] }  ~^T7^T'^L  * 

SinntjJq^  M-qjr 

as3  =  (M-q-2/  Bb  sinh(P2o)  tanhh((iq.2/)  -  [H20  Bb  cosh(^20)  +  -5] }  -C-Sh(^~2^  ^  , 

cosn^liq^  (Jqjr 

=  {*  Mq+2jr  Bb  sinh(H20)  tanh(pq+2j)  +  [M-20  Bb  coshto)  +  -5] }  COSh.^+2^  , 

cosntjiq^  |J^qj 

bs3  =  {Hq.2j  Ab  cosh((i20)  -  [M-20  Ab  sinh(^o)  -  -5]  tanh((iq-2j))  ~^^((^>r)  ^  , 

bs4  =  {-  Hq+2/  Ab  cosh(ji2o)  +  [M-20  Ab  sinh(^o)  -  -5]  tanh(Hq+2,r)}  ^  > 

ac  1  =  -as  1, 
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ac2  = 

-as2. 

bc\  = 

-  fol. 

bc2  = 

-to2. 

ac3  = 

-as3, 

ac4  = 

- as4 , 

bc3  = 

-  bs3. 

bc4  = 

-  bs4, 

as6  = 

ac6-‘ 

^.Kq 

2 

i|tanh(Mq,r)  - 

1  ) 
HA 

bs5  = 

bc5  = 

— ^jcothfpqj) 

- 1 ) 
HA 

and 

as5  = 

ac5  = 

bs6 

=  bc6  =  Kq( 

U0  +  %! 
,  °  2> 

B.2  Coefficients  For  Topographic  Instability 


Coefficients  for  (4.4.4)  -  (4.4.6)  axe  as  follows: 

al=XK(tanh(H20)-^ 

1 

M-20 


b2  =  A.K|coth(H20)  -  ^ 

-2KIU°  +  |), 


b\  =  a.2 


ul  =  K  (-7^—  cosh(|i20)  + 

\4|X20  2^20/, 


“2=K(-*S“+^1 

a3  =  .2Kf(coShto).^) 

b3  =  2K—  sinh(|i20) 

^  j 

and 

m3  =  0. 
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